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OPTIMIZATION  MODEL  FOR  LAND  TREATMENT  PLANNING,  DESIGN  AND  OPERATION 
PART  III.  MODEL  DESCRIPTION  AND  USER’S  GUIDE 


Jaclyn  A.  Baron  and  Daniel  R.  Lynch 

INTRODUCTION 

LTMOD  is  a  Fortran-coded  model  that  generates  optimal  combinations  of 
design  and  monthly  operating  parameters  for  slow-rate  land  treatment  of 
wastewater  systems  with  a  given  quantity  of  incoming  effluent.  It  is  a 
generalized  formulation  of  the  models  developed  in  Evaluation  of  Design  and 
Operating  Options  for  Slow-Rate  Land  Treatment  Systems  in  Cool/Humid  Areas 
(Baron  1981). 

LTMOD  consists  of  a  main  program,  which  coordinates  several  input  and 
output  files,  and  a  subroutine  containing  a  set  of  equations  representing 
the  physical  relationships,  the  groundwater  quality  constraints,  and  the 
objective  function  of  the  land  treatment  system.  These  equations  are 
coupled  with  additional  subroutines  comprising  a  generalized  reduced- 
gradient,  nonlinear  programming  algorithm.  User  input  is  required  for  1) 
parameters  of  the  LTMOD  equations — wastewater  influx,  physical  constants, 
desired  environmental  quality,  and  design  objective  for  the  site,  and  2) 
properties  of  the  optimization  procedure,  including  constraint  types, 
variable  bounds,  an  initial  solution  point,  and  various  tolerances  within 
the  algorithm.  Together,  these  two  input  types  define  the  physical  system, 
the  operating  options,  and  the  feasible  solution  space. 

This  report  defines  the  model  capabilities  and  underlying  assumptions, 
presents  an  overview  of  the  model  structure,  details  the  LTMOD  equations 
and  their  parameters,  discusses  the  procedure  for  running  the  model,  and 
includes  a  sample  problem  and  results.  Input  for  the  optimization 
procedure  is  described  only  when  it  is  important  that  it  be  consistent  with 
the  LTMOD  equations.  Additional  options  and  detailed  formatting  specifica¬ 
tions  for  input  for  the  optimization  procedure  are  in  the  GRG2  User's  Guide 
(Lasdon  et  al.  1978),  with  which  this  manual  should  be  used. 

MODEL  CAPABILITIES  AND  UNDERLYING  ASSUMPTIONS 

LTMOD  recognizes  the  dual  nature  of  the  land  treatment  process  whereby 
nitrogen  in  the  effluent  is  removed  by  biological,  physical  and  chemical 
processes  in  both  the  storage  facility  and  in  the  soil-crop  system.  The 
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Figure  l.  Land  treatment  design 
alternatives. 

resultant  design  tradeoff  between  storage  capacity  and  irrigation  area  is 
the  fundamental  issue  behind  the  model.  A  general  representation  of  the 
full  set  of  design  alternatives  is  shown  in  Figure  1.  The  feasible  design 
options  are  bounded  by  1)  AO,  the  minimum  area  to  which  effluent  can  be 
applied  to  avoid  ponding  (based  on  the  drainage  properties  of  the  site)  ,  2) 
SO,  the  minimum  storage  capacity  required  to  hold  effluent  during  the 
months  when  application  is  not  permitted,  and  3)  combinations  of  irrigation 
area  _>  AO  and  storage  capacity  SO  that  are  just  sufficient  to  meet  the 
groundwater  quality  constraints.  The  lowest  cost  designs  lie  along  or  at 
the  ends  of  segment  0,  at  the  edge  of  the  feasible  region.  Based  on  the 
potential  for  renovation  of  the  effluent  nitrogen  in  storage  and  on  the 
land  and  the  permissible  nitrogen  concentration  in  the  groundwater  under¬ 
lying  the  site,  LTMOD  may  be  used  to  satisfy  one  of  three  design 
objectives: 

1)  Find  the  minimum  required  effluent  storage  facility  capacity  S 
that  is  associated  with  a  given  irrigation  area  A  where  A  AO. 

2)  Find  the  minimum  required  irrigation  area  A  that  is  associated 
with  a  given  effluent  storage  facility  capacity  S  where  S  SO. 

3)  Find  the  minimum  required  irrigation  area  AO  and  the  minimum 
effluent  storage  facility  capacity  S^o  associated  with  this 
area. 

The  quantity  of  nitrogen  removed  from  the  effluent  in  storage  and  on 
the  land  (and  in  the  system  as  a  whole)  depends  on  dynamic  processes  and 
thus  on  the  operating  characteristics  of  the  site.  Operating  decision 
variables  include  the  volume  of  effluent  stored,  the  volume  of  effluent 
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applied  to  the  irrigation  area  from  storage,  and  the  volume  of  effluent 
that  bypasses  the  storage  facility  for  each  month  throughout  one  year  of 
operation.  Additional  system  operating  properties,  which  follow  from  these 
controllable  operating  decisions,  include  the  drainage  from  the  site,  the 
nitrogen  concentration  in  the  percolate,  the  nitrogen  concentration  in 
storage,  and  crop  nitrogen  uptake  for  each  month. 

The  LTMOD  equations  explicitly  or  implicitly  include  certain  assump¬ 
tions  that  delineate  the  type  and  properties  of  the  system  to  which  the 
model  may  be  applied.  Specifically,  it  must  be  a  single-crop  system  in  a 
cool  and/or  humid  area,  and  the  groundwater  nitrogen  concentration  must  be 
the  limiting  environmental  concern. 

The  model  assumes  that  a  single  crop  is  grown  on  the  entire  irrigated 
site,  that  the  cropped  area  is  equal  throughout  its  growing  season,  and 
that  the  entire  area  is  evenly  irrigated.  The  planting  and  harvesting 
regime  is  exogenously  fixed. 

The  model  is  for  a  cool  and/or  humid  area.  The  mass  balances  by  which 
the  drainage  from  the  site  is  calculated  are  independent  for  each  month  and 
do  not  provide  for  changes  in  the  soil  moisture  content.  It  is  assumed 
that  the  soil  is  always  kept  at  or  near  field  capacity.  It  follows  that 
sufficient  soil  moisture  is  available  to  the  crop  at  all  times  to  avoid 
water  deficit  and  stress.  Crop  moisture  uptake  is  assumed  to  be  inde¬ 
pendent  of  effluent  application  and  to  equal  the  potential  evapotranspira- 
tion  rate  in  each  month. 

Groundwater  nitrogen  concentration  is  assumed  to  be  the  limiting 
environmental  consideration.  Nitrogen  is  the  only  system  contaminant 
explicitly  modeled.  The  concentration  of  nitrogen  in  the  percolate  from 
the  site  may  be  constrained  on  an  annual  and/or  monthly  basis.  Nitrogen 
removal  in  the  storage  facility  is  represented  by  an  overall  monthly 
first-order  reaction  constant.  Monthly  crop  uptake  is  represented  by 
exponential  functions  of  the  nitrogen  available  to  the  crop.  The  monthly 
nitrogen  mass  balances  are  interdependent  and  allow  for  the  soil  storage 
and  release  of  applied  nitrogen  in  succeeding  months  that  has  been 
observed,  especially  during  winter  application  of  effluent  at  experimental 
sites.  This  is  accomplished  by  estimating  the  percentage  of  the  nitrogen 
in  the  soil  system  that  is  leached  in  each  month  (not  including  native  soil 
nitrogen,  which  is  assumed  to  remain  constant). 
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OVERVIEW  OF  MODEL  STRUCTURE 

The  model  consists  of  the  sections  shown  in  Table  1.  Figure  2  shows 
the  interactions  between  the  model  components.  The  main  program  allocates 
memory  and  numbers  the  input  and  output  files.  It  calls  subroutine  GRG , 
which  is  part  of  the  GRG2  file,  where  the  optimizing  process  occurs.  The 
INPUT  and  DEFINE  files  are  read  in  the  DATAIN  subroutine  and  GCOMP  is 
called  upon  as  needed.  The  PARAM  file  is  read  by  GCOMP  on  the  initial 
iteration.  The  output  is  sent  to  the  file  specified  for  this  purpose  by 
various  subroutines  in  the  GRG2  program. 

The  LTMOD ,  GRG2  and  DEFINE  files  must  be  available  and  accessible  on 
the  data  processing  system.  The  GRG2  file,  the  DEFINE  file  and  the  GCOMP 
portion  of  the  LTMOD  file  are  used  as  they  are;  the  MAIN  program  may 
require  modifications  pertaining  to  the  allocated  memory  and  file 
numbering,  depending  on  the  type  of  system  for  which  the  model  is  used. 

Table  1.  LTMOD  files. 


File  Name 


Contents 


LTMOD 

GRG2 

DEFINE 

PARAM 

INPUT 

OUTPUT 


Main  program  and  subroutine  GCOMP  (land  treatment  equations). 
Nonlinear  optimization  algorithm  subroutines. 

Names  of  LTMOD  functions  and  variables  for  labeling  output. 
Parameters  of  LTMOD  equations. 

Variable  and  function  types,  variable  bounds  and  options  for 
optimization  procedure. 

Output  file  of  intermediate  information  and  results. 


LTMOD  GRG2 


Figure  2.  Model  components  and  their  interactions. 
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The  INPUT  and  PARAM  files  are  supplied  by  the  user.  All  files  may  be 
arbitrarily  named. 

LTMOD  FILE 

The  generalized  land  treatment  formulation  contained  in  subroutine 
GCOMP  in  the  LTMOD  file  consists  of  85  functions  comprising  the  system 
constraints,  an  objective  function,  and  7  calculations  of  auxiliary  system 
properties  that  are  not  Included  in  the  optimization  process.  The  func¬ 
tions  include  85  decision  variables  and  160  parameters.  GCOMP  need  not  be 
modified  by  the  user;  however,  an  understanding  of  the  land  treatment 
equations  is  necessary  to  develop  the  user-supplied  PARAM  and  INPUT  files 
and  to  interpret  the  model  results.  A  schematic  of  the  system,  including 
the  major  variables  and  parameters,  is  shown  in  Figure  3.  Below,  the  93 
individual  equations  are  described,  using  the  equation  numbering  scheme 
G ( 1 )  through  G(93)  that  is  used  in  the  GCOMP  code. 

Note  that  continuous  multiyear  operation  of  the  system  is  assumed,  and 
thus  all  variables  defined  at  the  end  of  month  12  are  required  to  be 
identical  to  those  at  the  beginning  of  month  1.  Any  occurrence  of  the 
subscript  0  implies  the  subscript  12. 

Nitrogen  Crop  nitrogen 

renovation  uptake 


Effluent  bypass  Drainage 

Figure  3.  Simplified  schematic  of  a  land  treatment 
system.  Variables  are  defined  in  the  text. 

LTMOD  Constraints 

Equations  G( 1 )  -  G( 12)  represent  the  monthly  mass  balance  of  water  in 
storage: 

0 ,  +  B,  +  S,  -  S,  ,  +  (E,  -  P,  )  »max( S  ...S,_)/DP  =  I,  (1) 

kkkk-1  kk  1  12  k 
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where  the  decision  variables  are 

»  volume  of  effluent  applied  from  storage  In  month  k 
Bk  *  volume  of  effluent  applied  which  bypasses  storage  in  month  k 
Sk  *  storage  volume  at  end  of  month  k 
and  the  parameters  are 

Ek  »  evaporation  volume  in  month  k  per  unit  area 

Pk  *  precipitation  volume  in  month  k  per  unit  area 

DP  =*  depth  of  storage  facility  (design  depth  at  maximum  volume) 

Ik  *  incoming  effluent  volume  in  month  k. 

Equations  G ( 1 3 )  -  G( 24)  represent  the  monthly  mass  balance  of  water  at  the 
irrigation  site: 

EFFk*(°k  +  V  +  (Pk  “  Wk),A  “  Dk  =  °  (2) 

where  the  decision  variables  are 

(\  =  volume  of  effluent  application  from  storage  in  month  k 
Bk  =  volume  of  effluent  application  which  bypasses  storage  in  month 
k 

A  *  irrigation  area 

=  volume  of  percolate  from  site  in  month  k 
and  the  parameters  are 

EFFk  =  irrigation  efficiency  in  month  k  including  runoff 
Pk  =  precipitation  depth  in  month  k 
Wk  =  crop  water  uptake  in  month  k  (depth) . 

Equations  G(25)  -  G(36)  require  that  the  soil  drainage  capacity  in  each 
month  must  not  be  exceeded: 

MDR.  .A  -  D,  >  0  (3) 

k  k  — 

where  the  decision  variables  are 
A  ■  irrigation  area 

»  volume  of  percolate  from  site  in  month  k 
and  the  parameter  is 

MDRk  ■  soil  drainage  capacity  in  month  k  (depth). 

Equations  G(37)  -  G(48)  represent  the  monthly  mass  balance  of  nitrogen  at 
the  irrigation  site: 

l/FNL^*(NDk)  -  (1  -  DeNk).(Ok.NCAk  +  +  (1  -  FNL^^/FNL •NDk_1 1 

+  A«NUk  -  0  (4) 
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where  NCA^ ,  the  average  nitrogen  concentration  in  storage  in  month  k,  is 
approximated  as 


NC\-<%-\  +  NCk-r\-i)/<sk  +  sk-i) 


(5) 


The  decision  variables  are 

ND^  =  mass  of  nitrogen  leached  to  groundwater  to  month  k 
0^  =  volume  of  effluent  application  from  storage  in  month  k 
NC^  =  nitrogen  concentration  in  storage  at  end  of  month  k 
=  storage  volume  at  end  of  month  k 

=  volume  of  effluent  application  which  bypasses  storage  in  month 
k 

A  =  irrigation  area 

NU^  =  crop  nitrogen  uptake  in  month  k  (mass/area) 
and  the  parameters  are 

FNL^  =  fraction  of  soil  system  nitrogen  imbalance  leached  in  month  k: 
ND^  =  [  FNL^/ ( 1 -FNL^)  ]  •  NSSj^  ,  where  NSS^  is  the  nitrogen 
stored  in  the  soil  system  at  the  end  of  month  k 
(1-FNLfc)  =  fraction  of  soil  system  nitrogen  imbalance  retained  in  soil 
storage  at  the  end  of  month  k 
DeN^  =  fraction  of  nitrogen  in  soil  system  lost  through 
denitrification  and  volatilization  in  month  k 
NI^  =  nitrogen  concentration  of  incoming  effluent  in  month  k. 
Equations  G(49)  -  G(60)  describe  the  monthly  crop  uptake  of  nitrogen, 
represented  as  an  exponential  function  of  the  nitrogen  applied.  Its  basic 
form  is 

NU,  =  CP,  •  MNU,  .[l  -  EXP(-b/MNU,  )]  (6) 

k  k  k  k 


where  b,  the  nitrogen  available  to  the  crop  per  unit  land  area  is 

b  =  [ N 1^ • Bk  +  NCAk-Ok  +  (1  -  FNLk_1)/FNLk_1 .NDR_1 ]/A  .  (7) 

Substitution  for  NCA^  and  rear rangraent  yields 

[(" 

,-i-w) 


NU,  -  CP,  »MNU,  .71  -  EXP 
k  k  k  I 


»VBk-U  -  ™Vl>/™W«Dk_1J.(Sk+  Sk_,) 


-  0,  .(NC,  «S,  +  NC, 
k  k  k  k- 


1  / 1 <  sk  +  Sk_! ) «A»MNUk ] 


=  0 


(8) 
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The  decision  variables  are 


NUk  =  crop  nitrogen  uptake  in  month  k  (raass/area) 

Bk  =  volume  of  effluent  application  which  bypasses  storage  in  month 
k 

NDfc  =  mass  of  nitrogen  leached  to  groundwater  in  month  k 
Sk  =  storage  volume  at  end  of  month  k 

0k  =  volume  of  effluent  application  from  storage  in  month  k 
NCk  =  nitrogen  concentration  in  storage  at  end  of  month  k 
A  =  irrigation  area  in  month  k 
and  parameters  are 

MNUk  =  maximum  crop  nitrogen  uptake  in  month  k  (mass/area) 

Nlk  =  nitrogen  concentration  of  incoming  effluent  in  month  k 
FNLk  =  fraction  of  soil  system  nitrogen  imbalance  leached  in  month  k 
CPk  =  1  if  crop  is  grown  in  month  k,  =  0  if  it  is  not.  CPk  must 

also  be  set  to  zero  if  crop  does  not  take  up  nitrogen  in  a  month 
of  growth. 

Equations  G(bl)  -  G(72)  represent  the  monthly  mass  balance  of  nitrogen  in 
storage: 


NC,  *S,  -  NC,  ,-S,  ,  +  0,  .NCA,  +  NK,  »NCA,  *SA, 

k  k  k~2  k-1  k  k  k  k  k 

+  NI,  •  B,  =1,  •  NI, 

k  k  k  k 

where  SA^ ,  the  average  storage  volume  in  month  k  is: 

SA,  =  (S,  +  S,  )  /2 
k  k  k-1 

Substituting  for  NCAk  and  SAk  and  simplifying  terms  yields: 

NC,  -S,  -  NC,  ,-S,  ,  +  0,  -(NC,  -S,  +  NC,  ,.S,  .  )/(S,  +  S,  ,  ) 
k  k  k-1  k-1  k  k  k  k-1  k-1  k  k-1 


+  2  "V^k^k +  NCk-rsk-i>  +  NIk*Bk  =  NIk-xk 


(9) 


(10) 


(11) 


The  decision  variables  are 

NCk  =  nitrogen  concentration  in  storage  at  end  of  month  k 
Sk  =*  volume  of  storage  at  end  of  month  k 

0^  =*  volume  of  effluent  application  from  storage  in  month  k 

Bk  *  volume  of  effluent  application  which  bypasses  storage  in  month 


and  the  parameters  are 

NK(^  =  percentage  of  nitrogen  removed  in  storage  in  month  k 
Nl^  =  nitrogen  concentration  of  incoming  effluent  in  month  k 
Ik  =  incoming  effluent  volume  in  month  k. 

Equation  G(73)  enforces  the  annual  percolate  nitrogen  concentration 
constraint . 

NCQ*£d  "  JND  >  0  .  (12) 

k  k  k  k 

The  decision  variables  are 

Dk  =  percolate  volume  from  site  in  month  k 
NDk  =  mass  of  nitrogen  leached  to  groundwater  in  month  k 
and  the  parameter  is 

NCQ  =  maximum  permissible  annual  average  nitrogen  concentration  in 
the  groundwater. 

Equations  G(74)  -  G(85)  enforce  the  monthly  percolate  nitrogen 
concentration  constraint: 

NCQM •  D,  -  ND,  >  0  .  (13) 

k  k  — 

The  decision  variables  are 

Uk  =  percolate  volume  from  the  site  in  month  k 
NDk  =  mass  of  nitrogen  leached  to  groundwater  in  month  k 
and  the  parameter  is 

NCQM  =  maximum  permissible  monthly  nitrogen  concentration  in  the 
groundwater  in  month  k. 

LTMOD  objective  function 

Equation  G(86)  states  the  objective: 

Minimize  [OA* A  +  OB »MAX( S!  ...S12)]  (14) 

where  the  decision  variables  are 
A  <*  irrigation  area 

Sk  =  storage  volume  at  the  end  of  month  k. 

Varying  the  values  of  the  parameters  OA  and  OB  enables  a  choice  of  three 
objectives: 

1)  To  minimize  the  irrigation  area  with  a  given  storage  capacity, 
set  OA  *  1  and  OB  *=  0, 

2)  To  minimize  the  storage  capacity  with  a  given  irrigation  area, 
set  OA  =  0  and  OB  =1, 
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3)  To  minimize  the  irrigation  area  based  on  hydraulic  properties  of 
the  site  and  to  find  the  minimum  required  storage  associated  with 
this  area,  set  OA  =  1  and  OB  =  fractional  value  ensuring  that  the 
storage  term  in  the  objective  function  is  less  than  the  area  term 
at  all  times  in  the  optimization  procedure. 

It  is  essential  that  the  classification  and  bounds  on  the  storage  and 
area  decision  variables  in  the  INPUT  file  (see  below)  be  consistent  with 
the  choice  of  objective. 

LTMOD  additional  calculations 

Equations  G(87)  -  G(93)  represent  auxiliary  calculations  concerning 
the  system  size  and  the  various  renovation  processes.  The  mass  of  nitrogen 
removed  in  the  storage  facility  is  calculated  in  G ( 8 7 ) : 

C(87)  =  Ynk  •  (NC  *S,  +  NC,  , -S,  ,)/ 2  (15) 

:  k  k  k  k— 1  k— 1 

k 

while  the  percentage  of  nitrogen  removed  in  storage  is  calculated  in  G(88)  : 

G(88)  =  G(87 )/  y_(NI  •  l  ) .  (16) 

k 

Equations  G(89)  and  G(90)  calculate  the  comparable  mass-  and  percentage- 
removal  figures  for  the  crop: 


G( 89)  =  A.^NU 
k  k 

(17) 

G(90)  =  G(89)/)'(NIk.Ik). 
k 

(18) 

Equation  G(91)  calculates  the  total  mass  of  nitrogen  removed  by  the  tvro 
processes : 

G(9 1 )  =  G(87)  +  G(88)  (19) 

while  G(92)  gives  the  overall  nitrogen  removal  percentage: 

G(92)  =  G(91)/£(NI  -I  ).  (20) 

k 

Finally,  the  maximum  weekly  wastewater  application  rate  is  computed  in  eq 
G(93)  : 

G(93)  -  MAX(01  ...0  )/(4.3  A). 


(21) 


PARAM  FILE 

The  PARAM  fiLe  contains  the  parameters  for  the  LTMOD  land  treatment 
equations  and  is  supplied  and  named  by  the  user.  A  total  of  160  parameter 
values  must  be  specified.  The  file  must  consist  of  26  lines  with  6  values 
on  each  and  a  27th  line  with  4  values.  The  values  are  input  in  fields  of 
10  and  are  formatted  F10.6.  They  are  inserted  in  the  order  shown  in  Table 
2. 


Table  2.  Order  of  input  for  parameter  values.  The  GRG2  User's  Guide 
(Lasdon  et  al.  1978)  states  that  best  results  are  obtained  when  the 
absolute  values  of  the  problem  functions  are  less  than  102,  so  parameters 
should  be  scaled  accordingly. 


Lines  Parameter 


1 

and 

2 

1(12) 

Incoming  effluent  in  each  month  (volume). 

3 

and 

4 

P(  12) 

Precipitation  in  each  month  (depth). 

5 

and 

6 

E(  12) 

Evaporation  in  each  month  (depth). 

7 

and 

8 

W(  12) 

Crop  water  uptake  (PET)  in  each  month  (depth). 

9 

and 

10 

MRD( 12) 

Soil  drainage  capacity  in  each  month  (depth). 

11 

and 

12 

NIC  12) 

Nitrogen  concentration  of  incoming  effluent  in  each 
month  (mass/volume). 

13 

and 

14 

NK( 12) 

Percentage  of  nitrogen  removed  in  storage  in  each 
month. 

13 

and 

16 

DeN( 12) 

Fraction  of  nitrogen  in  soil  lost  through 
denitrification  and  volatilization  in  each  month. 

17 

and 

18 

CP( 12) 

Growth/no  growth  indicator  for  crop  in  each  month  (0.0 
or  1.0). 

19 

and 

20 

MNU( 12) 

*Maximum  potential  crop  nitrogen  uptake  in  each  month 
(mass/ area) . 

21 

and 

22 

FNL( 12) 

♦Fraction  of  nitrogen  in  soil  leached  in  each  month. 

23 

and 

24 

NCQM( 12) 

Maximum  permissible  nitrogen  concentration  in  percolate 
in  each  month  (mass/volume) . 

25 

and 

26 

EFF( 12) 

Irrigation  efficiency  in  each  month. 

27 

NCQ 

Maximum  permissible  annual  average  nitrogen 
concentration  In  groundwater  (mass/ volume) . 

27 

OA 

Objective  function  area  term  coefficient  (0.0  or  1.0). 

27 

OB 

Objective  function  storage  terra  coefficient,  0.0  <  OB  < 

1.0. 

Depth  of  storage  facility. 

27 

DP 

*If  the  crop  is  not  grown  in  month  k,  MNU^  should  be  set  to  some  arbitrary 
positive  value  to  avoid  division  by  zero. 

+The  fraction  of  nitrogen  leached  in  each  month  must  be  set  to  a  positive  value 
to  avoid  division  by  zero. 


11 


INPUT  FILE 

The  input  file  is  supplied  and  named  by  the  user.  The  file  must  be 
consistent  with  the  instructions  and  format  specifications  in  the  GRG2 
User ' s  Guide ,  Section  5.  The  required  parts  of  this  file  are: 

1)  Line  1  (Number  of  variables  and  functions). 

2)  ROWS  section  (Function  type  and  right-hand  sides). 

3)  BOUNDS  section  (Variable  type  and  bounds). 

4)  INITIAL  section  (Initial  values  of  decision  variables). 

5)  FUNCTION  and  VARIABLE  section  titles.  (Title  lines  only,  in  this 
order.  The  remainder  of  these  sections  is  read  from  the  DEFINE 
file.) 

The  user  may  wish  to  reduce  or  expand  the  level  of  printout,  reduce 
computation  time,  or  improve  performance  by  activating  additional  options 
discussed  in  the  GRG2  User1 s  Guide ,  Section  5. 

The  required  sections  are  discussed  in  detail  below.  A  portion  of  the 
data  will  be  identical  for  any  run  using  the  LTMOD  model.  The  user  should 
note  which  inputs  are  fixed  (Indicated  with  an  asterisk.)  and  which  depend 
on  the  particular  application. 

Line  1 

The  number  of  variables  is  85*,  and  the  number  of  functions  is  93*. 
Neither  of  these  numbers  may  be  modified  by  the  user. 

ROWS  section 

For  each  of  eq  G(l)  -  G(93)  the  equation  type  must  be  specified 
according  to  the  following  letter  code: 

E  indicates  an  equality  constraint, 

G  indicates  a  "greater  than  or  equal  to"  constraint, 

0  indicates  the  objective  function, 

N  indicates  a  constraint  that  is  ignored  in  the  optimization. 

TabLe  3  indicates  the  proper  choice  of  type  code  as  well  as  the  proper 
value  for  the  right  side  of  each  equation.  Only  eq  G(73)-G(85),  which  deal 
with  the  form  of  the  environmental  constraints,  require  user  decisions  with 
respect  to  equation  type.  Equation  G(73)  is  set  greater  than  or  equal  to 
zero  (G)  if  the  nitrogen  concentration  in  the  percolate  is  to  be 
constrained  on  an  average  annual  basis,  and  ignored  in  the  optimization 
procedure  (N)  if  it  is  not.  Likewise,  eqs  G(74)  -  G(85)  are  set  greater 
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Table  3.  Function  types  and  right  sides  for  ROWS  section  of  Input  file. 


Function 

type 

Right  si< 

G(l)-G(12) 

E* 

*k+ 

G(13)-G(24) 

E* 

0.0* 

G(25)-G(36) 

G* 

0.0* 

G(37)-G(48) 

E* 

0.0* 

G(49)-G(60) 

E* 

0.0* 

G(6  1)-G(72) 

E* 

Nlk*Ik 

G(7  3) 

G  or  N 

0.0* 

G( 74 )-G( 85) 

G  or  N 

0.0* 

G(86) 

0* 

-  * 

G(87)-G(93) 

N* 

-  * 

*These  inputs  are  fixed. 

+i(c  =  incoming  effluent  in  month  k  (volume), 

NI^  =  nitrogen  concentration  of  incoming  effluent  in  month  k  (mass/ volume) . 

than  or  equal  to  zero  (G)  if  the  nitrogen  concentration  in  the  percolate  is 
to  be  constrained  on  a  monthly  basis,  and  ignored  in  the  optimization 
procedure  (N)  if  it  is  not. 

BOUNDS  section 

The  variable  types  and  bounds  depend  heavily  on  the  particular 
application.  The  assignment  of  variable  bounds  describes  physical 
constants  and  constraints  of  the  land  treatment  system  that  are  not 
addressed  by  the  system  parameters  and  defines  the  operating  options  that 
may  be  considered  in  the  search  for  an  optimal  solution.  It  is  important 
that  variable  bounds  be  consistent  with  the  system  parameters.  This 
section  includes  definitions  of  the  85  LTMOD  decision  variables  and  some 
guidance  for  assigning  types  and  bounds.  As  in  the  ROWS  section,  a  letter 
code  signifies  the  type  of  bounds  on  each  variable: 

E  indicates  that  the  variable  is  fixed  at  a  given  value, 

G  indicates  that  the  variable  must  be  greater  than  or  equal  to  a 

given  value, 

R  indicates  that  the  variable  is  restricted  between  two  given 
values. 

X(l):  Irrigation  area  A.  If  the  objective  is  to  minimize  the  irri¬ 
gation  area  with  a  given  storage  capacity  (OA  -  1,  OB  -  0)  or  to  minimize 
the  irrigation  area  based  on  hydraulic  properties  of  the  site  and  find  the 
minimum  required  storage  associated  with  this  area  (OA  -  1 ,  OB  -  fractional 
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value),  then  set  X( 1 )  greater  than  or  equal  to  (G)  zero.  If  the  objective 
Is  to  minimize  the  storage  capacity  with  a  given  Irrigation  area  (OA  =  0, 

OB  =1)  then  set  X( t )  as  a  fixed  variable  (E)  at  the  area  of  interest. 

X(2)-X(13):  Storage  volume  at  the  end  of  each  month  S^.  The 

storage  volume  at  the  end  of  each  month  must  be  higher  than  some  small 
positive  value  to  avoid  division  by  zero.  if  the  objective  is  to  minimize 
the  irrigation  area  with  a  given  storage  capacity  (OA  =  1,  OB  =  0),  then 
set  X( 2 ) — X( 13)  as  restricted  variables  (R)  between  the  small  positive  value 
and  the  storage  capacity  of  Interest.  The  model  solution  will  always 
contain  at  least  one  storage  volume  that  is  equal  to  this  storage  capacity. 

If  the  objective  is  to  minimize  the  storage  capacity  witli  a  given 
irrigation  area  (OA  =  0,  OB  =  1)  or  to  minimize  the  irrigation  area  based 
on  hydraulic  properties  of  the  site  and  find  the  minimum  required  storage 
associated  with  this  area  (OA  =  1,  OB  =  fractional  value),  then  set 
X(2)-X(13)  as  greater  than  or  equal  to  (G)  variables  with  a  lower  bound  of 
the  small  positive  value. 

X(14)-X(25):  Monthly  effluent  application  from  storage  (\.  Monthly 

effluent  applications  should  be  set  to  be  greater  than  or  equal  to  (G)  zero 
in  the  months  when  application  is  permitted  when  irrigation  capacity  is  not 
a  concern,  and  as  restricted  variables  (R)  between  zero  and  the  irrigation 
capacity  if  appropriate.  Effluent  application  is  permissible  in  months 
when  crops  are  not  grown.  In  months  when  effluent  application  is  not 
permitted  due  to  climatic  considerations,  the  effluent  application 
variables  should  be  fixed  (E)  at  zero. 

X(26)-X(37):  Monthly  percolate  from  site  D^.  Monthly  percolate 

variables  should  be  set  greater  than  or  equal  to  (G)  zero  in  all  months. 

X(38)-X(49):  Nitrogen  concentration  in  storage  at  end  of  each  month 
NC^.  Nitrogen  concentration  in  storage  should  be  set  greater  than  or 
equal  to  (G)  zero  in  all  months. 

X(50)-X(61):  Monthly  crop  uptake  of  nitrogen  NU^.  Crop  uptake  of 
nitrogen  can  be  set  greater  than  or  equal  to  (G)  zero  in  all  months.  In 
months  when  crops  are  not  grown  in  nitrogen,  uptake  is  fixed  at  zero  by  the 
CP^  parameter. 

X(62)-X(73):  Monthly  quantity  of  nitrogen  leached  to  the  groundwater 
NDj<.  The  quantity  of  nitrogen  leached  variables  should  be  set  greater 
than  or  equal  to  (G)  zero  in  all  months. 
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X(74)-X(85):  Monthly  effluent  application  which  bypasses  storage 
Bfc.  If  bypassing  the  storage  facility  Is  an  operating  option,  the  bypass 
variables  should  be  restricted  (R)  between  a  lower  bound  of  zero  and  an 
upper  bound  of  1^,  the  Incoming  effluent  volume  In  each  month  when  appli¬ 
cation  is  permitted.  If  bypassing  storage  is  not  a  viable  option,  the 
bypass  variables  should  be  fixed  (E)  at  zero. 

INITIAL  section 

Initial  values  of  the  decision  variables  must  be  chosen  for  each  model 
run.  The  initial  solution  point  need  not  be  feasible;  however.  If  the 
Initial  point  is  poor,  a  feasible  solution  may  not  be  reached  and  either 
the  problem  will  be  declared  infeasible  or  various  error  messages  (size, 
division  by  zero)  will  be  generated.  LTMOD  is  particularly  sensitive  to 
the  relative  values  of  the  irrigation  area  and  the  effluent  storage 
capacity;  it  is  relatively  insensitive  to  the  initialization  of  the  other 
decision  variables  as  long  as  they  are  within  reason.  If  the  program  does 
not  complete  the  run,  modify  the  initial  irrigation  area  X(l)  when 
minimizing  the  area  with  a  given  storage  capacity,  or  modify  the  initial 
storage  volumes  X(2)-X(13)  when  minimizing  the  storage  capacity  with  a 
given  irrigation  area,  and  try  again. 

Depending  on  the  initial  solution,  the  global  optimum  may  not  be 
reached  for  reasons  that  are  discussed  in  the  GRG2  User’s  Guide  (Lasdon  et 
al.  1978).  If  the  final  irrigation  area  or  maximum  storage  volume  Is  less 
than  the  Initial  point,  it  is  a  good  idea  to  check  that  it  is  indeed  the 
minimum  requirement  by  lowering  the  initial  estimate  to  below  this  value 
and  rerunning.  Another  possible  check  that  the  solution  is  optimal  is  by 
comparison;  run  several  cases  around  the  point  of  interest  and  check  that 
none  of  the  solutions  is  out  of  line. 

A  reasonable  initialization  procedure  is: 

1)  Distribute  the  effluent  application  X(14)-X(25)  evenly  over  all 
the  months  in  which  it  is  feasible  and  set  it  to  zero  for  the 
remaining  months. 

2)  Set  the  drainage  volume  in  each  month  X(26)-X(37)  equal  to  the 
effluent  application. 

3)  Calculate  the  nitrogen  leached  in  each  month  X(62)-X(73)  by 
multiplying  the  drainage  volume  by  the  maximum  permissible 
nitrogen  concentration  in  the  percolate  (NCQ  or  NCQM) . 
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4)  Set  the  nitrogen  concentration  in  storage  in  each  month 
X(2)-X(13)  as  half  the  incoming  effluent  value. 

5)  Estimate  the  crop  nitrogen  uptake  in  each  month  in  which  it  is 
grown  by  multiplying  the  nitrogen  concentration  in  storage  by  the 
effluent  application  and  using  the  crop  nitrogen  uptake 
exponential  functions. 

6)  Bypass  variables  X(74)-X(85)  tend,  in  final  solution,  to  be  at 
their  lower  bound  (zero)  or  their  upper  bound  (the  incoming 
effluent  volume),  depending  upon  the  application;  it  is  worth¬ 
while  to  rerun  the  model  with  low  and  high  initial  bypasses  to 
ensure  that  the  optimal  solution  is  found. 

7)  When  minimizing  the  irrigation  area  associated  with  a  given 
storage  capacity,  set  the  storage  in  each  month  X(2)-X(13)  at 
approximately  half  the  maximum.  Roughly  estimate  the  hydrau¬ 
lically  feasible  area  by  dividing  the  effluent  volume  applied 
each  month  by  an  average  drainage  capacity.  Set  the  initial  area 
X ( 1 )  somewhere  above  this  result. 

When  minimizing  the  storage  capacity  associated  with  a  given  irriga¬ 
tion  area,  find  the  storage  capacity  required  to  store  the  incoming 
effluent  volume  during  the  months  when  irrigation  is  not  feasible;  set  the 
storage  volumes  X(2)-X(13)  somewhere  above  half  of  this  result. 

RUNNING  LTMOD  ON  THE  PRIME  SYSTEM 

Running  LTMOD  as  a  segmented  run  file  on  a  Prime  system  is  straight¬ 
forward.  The  LTMOD  file  may  be  used  with  no  modification.  The  allocated 
memory,  PARAMETER  MAXMEM,  is  20000.  The  file  numbering  is: 

INPUT  LOGICAL  UNIT  5 

PARAM  LOGICAL  UNIT  7 

DEFINE  LOGICAL  UNIT  8 

OUTPUT  LOGICAL  UNIT  6. 

The  run  procedure  is  as  follows: 

1)  Compile  LTMOD  and  GRG2  files.  Either  the  Fortran  IV  (FTN)  or 
Fortran  77  (F77)  compiler  may  be  used: 

F77  LTMOD  or  FTN  LTMOD  -64V 

F77  GRG2  -BIG  FTN  GRG2  -64V  -BIG. 
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2)  Load  Segmented  Runflle: 

SEG  (initialize  SEG  utility) 

lo  #LTM  (load  runflle  name) 

lo  b-LTMOD  (load  compiled  LTMOD) 

10  b-GRG2  (load  compiled  GRG2) 

11  vapplb  (v-mode  applications  library) 

li  (standard  Fortran  library) 

save  (save  runflle) 

Q  (quit  SEG  utility  and  return  to  Primos  command  level). 

The  package  is  now  saved  for  all  future  runs.  To  run  in  foreground,  type 
SEG  //LTM  (or  runflle  name)  and  input  the  INPUT,  PARAM,  DEFINE  and  OUTPUT 
file  names  as  requested. 

A  typical  background  job  might  look  like: 

JOB  user  It  -HOME  'user  it  password'  -QUEUE  LONG 

COMO  RUNLIST 

SEG  It  LTM 

INPUT 

PARAM 

DEFINE 

OUTPUT 

LOGOUT 

Runtime  error  messages  are  available  by  listing  the  RUNLIST  file. 

SAMPLE  PROBLEM  AND  RESULTS 
Problem  statement 

The  sample  system  is  a  hypothetical  10-mgd  facility  with  climatic  and 
crop  data  typical  of  central  New  Hampshire.  The  incoming  effluent  volume 
is  evenly  distributed  over  the  year  (1,151,300  m'Vmonth)  .  Monthly  esti¬ 
mates  of  precipitation  and  evaporation  at  the  site  are  shown  in  Table  4.  A 
mixed  forage  crop  is  grown  on  the  site  from  early  April  to  mid-November, 
and  its  potential  evaporation  rate  is  assumed  to  equal  the  evaporation  rate 
throughout  its  growing  season.  The  soil  is  a  sandy  loam  with  a  maximum 
hydraulic  loading  of  8  in. /week  (87  cm/month).  Irrigation  efficiency  is 
70%  and  runoff  is  assumed  to  be  negligible. 

The  incoming  effluent  contains  40  mg/L  nitrogen.  The  nitrogen  removal 
<n  the  storage  facility  is  20%  in  each  month  from  April  through  November 
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Table  4.  Climatic  data  for  the  sample  problem. 


Month 

Precipitat ion 
(cm) 

Potential 

Evaporation 

(cm) 

January 

6.9 

0.0 

February 

5.8 

0.0 

March 

6.6 

0.0 

April 

6.6 

3.0 

May 

7.9 

7.7 

June 

8.6 

11.2 

July 

8.9 

13.4 

August 

8.9 

11.3 

September 

8.1 

7.6 

October 

7.9 

3.9 

November 

7.1 

0.6 

December 

6.6 

0.0 

Table  5.  Estimated  monthly  percentages  of  sea¬ 
sonal  nitrogen  uptake  for  the  sample  problem. 


Month _ %  Seasonal  N  Uptake 


April 

15 

May 

20 

June 

15 

July 

20 

August 

20 

September 

10 

October 

- 

November 

- 

and  10%  in  the  cold  months  from  December  through  March.  The  maximum 
seasonal  uptake  of  the  forage  is  470  kg/ha  under  a  harvesting  regime  that 
includes  three  cuttings.  The  percentage  of  the  total  seasonal  uptake  that 
may  be  expected  in  each  month  with  this  management  schedule  is  shown  in 
Table  5.  Maximum  monthly  nitrogen  uptake  values  are  estimated  by  multi¬ 
plying  the  maximum  seasonal  uptake  by  the  appropriate  percentages.  Nitro¬ 
gen  losses  by  denitrification  and  volatilization  are  assumed  to  be  negligi¬ 
ble.  All  nitrogen  applied  to  the  irrigated  area  in  the  effluent  that  is 
not  consumed  by  the  crop  is  assumed  to  leach  in  the  same  month  in  which  it 
is  applied. 

The  study  objective  is  to  find  the  minimum  required  irrigation  area 
associated  with  a  storage  facility  with  a  5,500,000-m3  capacity  and  12-ft 
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(3.66-m)  depth  to  achieve  an  annual  average  concentration  of  nitrogen  in 
the  percolate  of  10  mg/L.  Effluent  Irrigation  Is  feasible  only  during  the 
growing  season  of  the  forage  crop.  Bypassing  the  storage  facility  is  not 
permitted. 

Input  and  output  files 

The  PARAM  file  for  the  problem  is  shown  in  Table  Al .  Two  items 
illustrate  the  flexibility  available  to  the  user: 

1)  A  percentage  of  the  winter  precipitation  is  snow,  which  melts  and 
percolates  through  the  soil  in  the  spring.  In  this  example,  for  simpli¬ 
city,  100%  of  the  winter  precipitation  is  deferred  until  April  when  it  is 
assumed  to  enter  the  soil  system. 

2)  To  model  the  fact  that  effluent  application  is  permissible  only  in 
the  first  half  of  November,  the  drainage  capacity  is  halved  in  this  month. 

The  INPUT  file  for  the  run  is  shown  in  Table  A2.  The  complete  OUTPUT 
file  is  shown  in  Table  A3.  The  required  irrigation  area  is  1280  hectares. 
Refer  to  the  GRG2  User's  Guide  for  further  details  of  the  output  informa¬ 
tion. 
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APPENDIX  A:  INPUT  AND  OUTPUT  FOR  THE  SAMPLE  PROBLEM 


Table  A1 . 


11  513 

11.  513 

11  513 

11  513 

11  513 

11  513 

0.  0 

0.  0 

0.  0 

0  039 

0.  009 

0  081 

0.  0 

0  0 

0.  0 

0.  134 

0  133 

0  076 

0  0 

0.  0 

0  0 

0  134 

0  113 

0  076 

0.  87 

0  07 

0,  07 

0  87 

0.  37 

0  a7 

0.  04 

0  04 

0  04 

0  04 

0  04 

0  04 

0  10 

0.  10 

0  10 

0.  30 

0.  20 

0  30 

0  0 

0  0 

0  0 

0  0 

0  0 

0  0 

0  0 

0  0 

0  0 

1  0 

1.  0 

1  0 

0  01 

0.  01 

0  01 

0  00940 

0.  00940 

0  004  7 

1  0 

1  0 

1.  0 

1 . 0 

1. 0 

1  0 

0.  015 

0.  015 

0  015 

0  015 

0  015 

0  015 

0.  70 

0.  70 

0.  70 

0.  70 

0.  70 

0.  70 

0  010 

1.  0 

0.  0 

PARAM  /'le. 


1  1  513 

11  513 

11.  513 

11.  513 

11.  5)3 

11  513 

0.  335 

0.  079 

0.  C35 

0.  079 

0.  071 

0.  0 

0.  03 

0  077 

0.  112 

0.  039 

0.  006 

0  0 

0.  03 

0.  077 

0.  1  12 

0.  039 

0.  006 

0.  0 

0.  07 

0  37 

0  97 

0.  87 

0/435 

0.  87 

0.  04 

0  04 

0  04 

0  04 

0  04 

0  04 

0  30 

0.  30 

0  30 

0.  30 

0.  20 

0.  10 

0.  0 

0.  0 

0.  0 

0.  0 

0  0 

0.  0 

1.  0 

1.  0 

1.  0 

c.  0 

0.  0 

0  0 

0.  00705 

0.  00940 

0.  00905 

C.  01 

0.  01 

0.  01 

1.  0 

1.  0 

1 . 0 

1.  0 

1 . 0 

1. 0 

0.  015 

0.  015 

0.  01  5 

0  015 

0.  01  5 

0.  01  5 

0.  70 

0.  70 

0.  70 

0.  70 

3.  6576 

0  70 

0.  70 
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PREVIOUS  PAGE 
IS  BLANK 


k 


Table  A2.  INPUT  file 


S5 

73 

ROWS 

E  1 

12 

11  513 

E  13 

24 

0  0 

0  25 

26 

0.  0 

E  37 

48 

0  0 

E  47 

60 

0.  0 

E  61 

72 

0.  4  c  05 

G  72 

0  0 

N  74 

35 

0  86 

N  87 

73 

END 

BOUNDS 

G  1 

1 

0.  0 

R  2 

13 

5297 

E  14 

16 

0.  0 

G  17 

24 

0  0 

E  25 

0.  0 

G  26 

37 

0  0 

G  33 

49 

0.  0 

G  50 

61 

0.  0 

G  62 

73 

0  0 

E  74 

95 

0  0 

END 

INITIAL 

SEP 

1 

50  0 

2 

13 

30.  0 

14 

16 

0  0 

17 

24 

13.  9 

25 

0  0 

26 

23 

0.  0 

27 

36 

16  7 

37 

0.  0 

33 

47 

.  024 

SO 

52 

0  0 

53 

.  005 

54 

.  0065 

55 

.  005 

56 

57 

.  0065 

53 

.  0033 

57 

61 

0  0 

62 

64 

0  0 

65 

72 

.  187 

73 

0.  0 

74 

35 

0  0 

END 

FUNCTIONS 

VARIABLES 

PRINT 

IPR  0 

END 

C-0 
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Tablt-  AJ.  Ol'TPUT  file 


1 

ONUilbFR  G.r-  VAhlAELTS  13  B 5 
Ot .'UMBER  Or  FUNCTIONS  1C  93 

OCR  AC  E  EE~FRVcD  TCP  UES5IAN  HAS  DIMENSION  05 
OLIMIT  ON  BINDING  CONSTRAINTS  13  92 

OACTUAL  LENGTH  Ur  Z  ARRAY  IS  19075 


HOWS 


E 

1 

12 

1  15130CS-01 

E 

13 

24 

0. OOOOOOE-Ol 

G 

25 

36 

0  0000002-01 

E 

37 

4B 

o  oooooc 3-:  i 

E 

49 

60 

0  OOGOOOE-O 1 

E 

61 

72 

4. 6050GGE-01 

G 

73 

73 

0.  OOOOOOE-Ol 

N 

74 

85 

0 

96 

N 

87 

93 

FIND 

BOUNDS 

G  110  OOOOOGE-Ct 
R  2  13  5  299000S-C1  5  5CC0COE+O1 

e  14  i6  o  ecooocs-o  1 

G  17  24  0  CCCOCCE-Cl 

E  25  25  0  OOOOOC 2-01 

G  26  37  0  GCCGC  OE  -0- 1 

G  39  49  0  0CC0GCE-01 

G  50  61  0  OOvCCOE-OJ 

G  62  73  0.  COCC'O'H-  ;  i 

E  74  35  O  COOCOCE'Ql 

END 

INITIAL 

SEP 


1 

1 

5. 

ooecceoE 

♦01 

2 

13 

3 

OOOOOOCt ♦0 1 

14 

16 

0 

0 OOOOOC t 

-01 

17 

24 

1 

3 9 GO 00 72 

♦01 

25 

25 

0 

CCCOCC 72 

-01 

26 

28 

0 

COCO O  '  -  0 

-01 

29 

36 

1. 

Cr* 0  0 0 0 C  £ 

<-01 

37 

37 

0 

OOOCCOCC 

-01 

33 

49 

o 

c 

4 OOOOOC 0 

-02 

50 

52 

0 

OGCOCOCE 

-01 

53 

53 

5 

GGCOOCC  E 

-03 

54 

54 

6. 

500000CS 

-03 

55 

55 

5 

OOGCCOCE 

-03 

56 

57 

6 

5CC0CG :e 

-03 

53 

53 

-3 

3000CCC  E 

-03 

59 

61 

0 

oooooc :e 

-01 

62 

64 

0 

OOOOOOOE 

-01 

65 

72 

1 

C-QOO^C 

-01 

73 

73 

0. 

OOOOOOOE 

-01 

74 

05 

0 

OOOOOOOE 

-01 

END 

FUNCTIONS 
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VAR  I AuL  E3 


PRINT 
IPR  0 

END 
GO 


OEPNFWT  =  1.000CE-04  EPINIT  =  I  COOOE-04  EP5TCP 

EPP1V  =  1 .  OOGOE  03  PI II EPS  =  G  0G0CGE-G1 

ONSTOP  =  3  ITLIM  -  10  LIHSFR  =  10000 

0 1  PR  =  0  PN4  =  0  PN5  =  0  PN6  =  0  PER  = 

OTANGENT  VECTORS  WILL  2E  USED  FOR  INITIAL  EST I I1ATES 


0 

or 


1  OOGOE -04 


DUMP  «=  0 

BASIC  VARIABLES 


OTHE  FINITE  DIFFERENCE  PARSH  USING  FORWARD  DIFFERENCES  WILL  BE  USED 


OOBJECTIVE  FUNCTION  WILL  PE  MINIMIZED 


ODFP  WILL  BE 

USED  IF  !* 

51 

JPERBASI 

CCS  <  =  B5 

PARAMETERS 

11. 513999 

11 

513999 

1 1 

.  513999 

1 1 .  513999 

11  513999 

1  1 

512999 

11  513979 

1 1 

51-999 

1 1 . 

512999 

1 1  513999 

11  513999 

1  1 

512999 

0  000000 

0 

000000 

0 

000000 

0.  325000 

0.  079000 

0 

csgooo 

0.  039000 

0. 

C5e000 

c* 

C  3 1 000 

0. 079000 

0.  071000 

0 

000000 

0.  000000 

0 

030000 

ocooco 

0  030000 

0.  077000 

0 

1 12000 

0  130000 

0 

123000 

0 

076000 

0  029000 

0.  006000 

0 

OCOOGO 

0.  ooocco 

0 

000000 

n 

0C0OQ0 

0  02 0000 

0  07/000 

0 

1 13000 

0.  130000 

0 

1 13000 

076000 

0  039000 

0.  006000 

0 

000000 

0.  0  70000 

0. 

S 70 COO 

c 

S70000 

0.  370000 

0.  370000 

0 

270000 

0.  070000 

0 

070000 

Q 

370000 

0  370000 

0.  435000 

0 

9 70000 

0.  000000 

0 

04C000 

0. 

040000 

0  04 0000 

0.  0-1 0000 

0 

040000 

0  040000 

0. 

040000 

r, 

040000 

0.  040000 

0.  04  0000 

0 

040000 

0  100000 

0 

100000 

0 

100000 

0  200000 

0  300000 

0 

200000 

0  300000 

0 

200000 

0 

300000 

0  300000 

0  200C00 

0 

1 00000 

0.  oococo 

0 

000000 

000000 

0.  000000 

0.  000000 

0 

000000 

0  000000 

0 

000OC0 

0 

000000 

0  oococo 

0  000000 

0 

00000  0 

0. 000000 

0. 

000000 

0 

000000 

1  000000 

1 .  000000 

1 

000000 

1 .  000000 

1 . 

000000 

1 

OOGOOO 

0.  000000 

0.  000000 

0 

000000 

0  010GC0 

0 

0 1 0000 

0. 

010000 

0  007050 

0  009400 

0 

007050 

0  009400 

0 

009400 

0 

004700 

0  010000 

0  0100:0 

0 

0 !  ('000 

1  000000 

1 

000000 

1 

000000 

1  .  000000 

1 .  OCOOGO 

1 

OCOOGO 

1.  000000 

1 

000000 

1. 

000000 

1 .  000000 

1  000000 

1 

000000 

0.  015000 

0. 

015000 

0 

015000 

0.  015000 

0.  015000 

0 

015000 

0.  015000 

0 

015000 

0 

015000 

0  oirooo 

0  01 5000 

0 

01 5000 

0.  700000 

0 

700000 

0. 

700000 

0  700 COO 

0  700000 

0 

700000 

0.  700000 

0 

700000 

0 

700000 

0  700000 

0.  700000 

0 

"COO 00 

0.  010000 

1 . 

000000 

0. 

000000 

3  657600 

1  OUTPUT  OF  INITIAL  VALUES 
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OSECTIGN  1  --  FUNCTIONS 

FUNCT  I  C'N 

INI T IAL 

LOWER 

UPPER 

NO  NAME  STATUS  TYPE 

VALUE 

L  I  N  J  T 

L  I M I  T 

1 

l./  0001 !  i 

5t  9 

C 

0 

■  A"  E  -  71 

1 

1 51 20GOC <01 

1  1 5130002-01 

t.AC.UL-iJ 

»  «  * 

r  g 

0 

:•••*•  f.  -M 

1 

15120002-01 

1.  151 3GCCC- 0 1 

3 

LA  GO  ON  3 

*  ^  ♦ 

f  Q 

0 

C'  OGOE  -01 

1 . 

1  513000E  ~<J1 

1  151 2CCCE  +  0 

1  *-4  V  »  _  7  . 

*<  » 

t  0 

1 

t-'C  022'  'E  '■0  1 

1 

151 2000S -C  1 

1  .  15130..  '"-01 

\j 

L*  CJl'.S 

<  *  * 

LG 

1 

1 

1 3 1 3G0QE -0 1 

1 .  1 513 oc  :r.  <-01 

•\ 

L  A' T  v 

f  * 

r  o 

1 

R  1  1 2255E +0 1 

1 

! 5 1 2000E -0 1 

1  151 300CE -0 1 

7 

LAs. .A."*  ' 

X  3  -a 

L>i 

1 

°2o9u''-'>C  <-01 

1 

! D 1 3C00E-0 1 

1  i5i3oc;r.  »oi 

Cl 

i  a  :  -s 

a  O  -4 

ro 

1 

c’260E,'#::t  *01 

1 

151  20001-  -01 

1.  1513C00E <01 

'  f 

LA  -L-.J-, 

Jt  A  » 

L  Q 

1 

ro5  :99  •<=  *■■:■  i 

1. 

1 0 1 3000E ♦O 1 

1.  i5i3o:.:: -01 

iO 

L  ACCOM  0 

a  ♦ 

r  q 

1 

3  7  *  n  ^  vf:  i 

1 

i 5 : 00002 -0 1 

1  1  51 3CC ( 2 -01 

1 1 

LAGOON  1  1 

•a  *  * 

r  o 

1 

-::.v r  l  *  a  -  o  i 

1 

161 JCOOE+O 1 

1 .  151 20c  :e  <0 

1.0 

LACCQN 1 2 

a  *•* 

EG 

0 

CUOC'CCOE  -  0  1 

1 

1 5 1 2000E -0 1 

1.  1513CCCE-0 

13 

H20  3L1 

EQ 

EO 

0 

CuC'OOOOE-O  1 

0 

GooooooE-oi 

0.  oooccoc <:-o 

14 

H20  DL.2 

EO 

EQ 

0 

OUOvOOO'C  -0  1 

0. 

CGOOOOOiZ— 0 1 

0  oooooccs -0 

IS 

H20  DL3 

EQ 

EQ 

0 

C  00000*  <E  -01 

0 

ouooocoE-oi 

0.  OOOOOC  .E-U 

16 

H20  EL 4 

>  *  * 

EQ 

9 

07V9R22E+00 

0. 

ooocoooE-oi 

0.  0000000*3-0 

17 

H20  CL 5 

<*•  * 

EO 

“5 

57000! »E»CO 

0 

OOOCOCOE-O 1 

0  OOOOOC  CE -0 

10 

H20  CL6 

*  41  * 

EO 

•‘6 

=7000 1 OE <00 

0. 

OOOCOOOE-Gl 

0.  OOOOCGCl-O 

17 

H20  GL7 

a  <>•» 

EG 

-7 

9.20000  /F  <00 

0 

fiCGCCOOE-Ol 

0.  oooocc-:.:-  c 

70 

H23  ELS 

EQ 

-ft 

570001  1.E  rC-0 

0 

0  'C-GOO>~E-0  1 

0  000000:1: 

01 

H2G  3L9 

EQ 

-5 

4200'.*:  IE  <-00 

0 

COG‘2GOCt2  -C  1 

0  ooooocoL 

•n  -•> 

CL  - 

H20  CLIO 

a  c  -a 

EG 

-3 

i*7  COu  1  *•£■*■  00 

0 

OOOCOOOc-O 1 

0  OOOOCCOC-O 

03 

H23  3l 1 1 

•5  L  <» 

ro 

_  'O 
c_ 

-i.OOC  k ;  ^  6C  *(.  0 

0 

OCCCCOOE-O 1 

0. OOOOCC C 3 

"T  /i 

H2J  ELI  2 

Cv< 

EO 

r. 

;ooc  **.  t:  -oi 

0 

COOCOOOE-O 1 

0. oooooc :c-g 

OS 

DR A  CPI 

GE 

4 

24997 44E+01 

0 

OOCCOGOE-Ol 

NONE 

26 

DR  A  CP  7: 

GE 

4 

349999  E-0 1 

0 

O-jC'C  OOOE-O) 

NONE 

07 

DRm  *.R3 

GE 

a 

3499794F+01 

0 

00  ?0(>0C,E -0  J 

NONE 

oc 

DR A  CP 4 

GE 

c.  . 

4  599994  E  *-0 1 

0 

0000 20 02-0 1 

NONE 

29 

DR A  CP 5 

GE 

~V 

C. 

45999‘>4E-»01 

0 

OOOCOOOE -0 1 

NONE 

20 

DR A  CP6 

GE 

—1 

C- 

.;5999-.:H"+01 

0. 

OGOOOOOE-Ol 

NONE 

31 

DRA  CF7 

GE 

n 

4 59974.  »E -01 

0. 

OCOOOGGE-Ol 

NONE 

32 

DRA  CP3 

GE 

T 
C.  . 

4  59 49 4  4  E-<- 01 

0. 

OOOOOOOE-O 1 

NONE 

13 

DRA  CP9 

GE 

45999? <•? ,o 1 

0 

Cr-~ OOOOH-'j  J 

NONE 

2  1 

DRA  CP  10 

•”'E 

2 

4  549994  E-OJ 

0 

orooocoE-ci 

NONE 

-3  r, 

DRA  CP  1 1 

GE 

<7> 

«A«V,?i?7  !  L+C  O 

0 

COCCOOCE-Ol 

NONE 

'-i  t 

DRA  CP  1 2 

r  l: 

4 

2  - 49944E  +  0 1 

0. 

OGGCOOOE-G 1 

NONE 

37 

N  DALI 

r  0 

f-G 

o<  c oe-o 1 

0 

GOOOOGOE-O 1 

0. OOOOCGC2-01 

30 

N  DAL 2 

EQ 

1  0 

0 

C  uC  OC’OOE-0 1 

0. 

OC  CCOOOE-O I 

0  oooocc: E- 0 1 

39 

N  DAL  3 

EQ 

1  0 

0 

COO  00002-0  1 

0 

CCCCCGOE-Q l 

0. OOOOCOC^-O 1 

40 

N  DAL 4 

r  o 

-I 

4  0  •  •  ■’  •  •<  f p 

0 

(>00 .  CCOE-O  1 

0  oooooocr  01 

41 

N  DAL5 

*  *r*- 

EG 

6 

0')0'.‘00(>E'OC 

0 

OC  's. OOOE-O  1 

0  OOOOOOOE-O 1 

42 

N  DAL 6 

it  ■* 

EQ 

-1 

4  6 0 C 0 0  0 E  -  0  <2 

0 

0-  'COOOE-Ot 

0  OOOOOOCC-Ol 

43 

N  DAL  7 

tt  ■* -fr 

EO 

ft 

O-IOOOCOE  -02 

0 

OOOCOOOE-Gl 

0.  ooooocc;:~gi 

44 

N  UAL 3 

■a  *  * 

EG 

ft 

OTOOOOOE  -02 

0 

00 1 OOGOE-O 1 

0.  OCOOOOOE-s)l 

46 

N  DAL 9 

a  i-  a 

EG 

-ft 

R60.0000C-02 

0 

OOOOOOGt-O 1 

0.  00000000-01 

46 

N  DAL  10 

*  »  ■* 

EG 

_  o 

o'*  ft 00 C  I'tZ  -0  l 

0 

Ou'OOCOE -0 1 

0  OOOOOC C2-0 1 

47 

N  DALI  1 

* 

t  0 

—  n 

CL 

446000(2-0  1 

0 

OOOOOOGE-O 1 

0  OOOOOOGE-O 1 

40 

N  DAL  12 

EG 

EQ 

0 

OOC  OOCOG-O  1 

0 

GCOOOOOE-Ol 

0  OOOOOOC  0-7*1 

47 

N  CROP  1 

EG 

EQ 

0 

C  000000*0  ~01 

0. 

COCCOOOE-O 1 

0  OOOOOOC 0-01 

50 

N  CROP 2 

EG 

EQ 

0 

OOOOCCA.E -01 

0 

OCOOOOCE-O 1 

0.  OOCOOOCO-Ol 

51 

N  CROP 3 

EG 

EQ 

0 

OOOOOOOE-O 1 

0. 

OOOGOOOE-O 1 

0  OOOOOCOO-Ol 

52 

N  CROP 4 

fr  r  * 

C  Q 

-1 

0312?2v'E  -04 

0 

OOOC-OOOE -01 

0  OOOOOOCO-Ol 

53 

N  CROP 5 

EQ 

6 

3036133E-04 

0 

OOOOOOOE-O 1 

0,  OOOOOC CO- 01 

54 

N  CR0P6 

*** 

EQ 

- 1 

03 1 372QE-04 

0 

OOOCOOOE-O 1 

0.  OOOOOOCE-O 1 

25 


or. 

N  CRGP7 

to 

6 

80B61332-04 

0  OOOOOOOE-O 1 

0  OOOOOOOE-O 1 

56 

N  CRCi-8 

tt  <  * 

LQ 

6 

8036 ! 72  C -04 

0. GQOCOOOE-G 1 

<)  0O0000CE  -01 

57 

N  CR0P9 

EQ 

-7 

1794029F.-O4 

0.  OOGCOOOE-0 1 

0  OOOCOOOE-Ol 

50 

N  CROP  10 

EQ 

EQ 

0. 

COOOOOOC -01 

0.  OGOOOOCE-Ol 

0  COOOOOOE-Ol 

09 

N  CROP  1  1 

BQ 

EQ 

0 

0000000£-01 

0. OOOOOOOE-O 1 

0.  OOOOOOOE  -O  1 

60 

N  CR0F12 

EQ 

EQ 

0. 

OOOOOOOE-O 1 

0.  OOOOOOOE-O 1 

0.  OOOOOOCE-O 1 

61 

LAG  Ml 

*  *?  •» 

EQ 

7 

199999oE-C2 

4  6050000E-G l 

4.  6050C0CC-0 ; 

62 

LAG  M2 

EQ 

7 

1  999996E-02 

4.  6050000E-0 1 

4.  605000CE-0 1 

63 

LAG  N3 

*  i.+ 

EQ 

7 

199999+E-02 

4.  6050000E-0 1 

4  6050000E -0 1 

64 

LAG  N4 

&  u  * 

EQ 

5 

9759999E-01 

4.  605C000E-0 1 

4, 60500CCE-0 1 

65 

LAG  NS 

EQ 

5 

9759999E-0 1 

4.  6050000E-01 

4  60500002  01 

66 

LAG  N6 

**.  * 

EQ 

5 

9759997E-0 1 

4  6GSCOOOE-0 1 

4.  6050GCC  2  -  0 1 

67 

LAG  N7 

EQ 

5 

9759999E-0 1 

4.  605CG00E-01 

4  60500:  ct;-oi 

63 

LAG  NO 

EQ 

5 

9759999E-0J 

4  605000CE-01 

4. 60500002-01 

69 

LAG  N9 

EQ 

5 

9759999E-01 

4  6G5C OOCE-0 1 

4  60  500001: -01 

70 

LAG  N 1 0 

EQ 

5. 

9759999E-01 

4  6050000E-0 1 

4.  6050000E-01 

71 

LAG  N 1 1 

*<  ■* 

EQ 

5. 

9759999E-01 

4.  60500C0E-0 1 

4  60500C  OE-Q 1 

72 

LAG  N 1  2 

<k  le* 

EQ 

7 

1 999996E-02 

4.  6050000E-01 

4  60500002 -0 1 

73 

GW  QUAL 

LL 

GE 

-1. 

4210055E -14 

0  ooooocoe-oi 

NONE 

74 

GWN1 

NA 

0 

OOOGOOOE -0 1 

75 

GW‘,'42 

NA 

0 

000000CE-01 

76 

GUN3 

NA 

0 

OOOGOOOE -01 

77 

GUN  4 

NA 

9. 

4499976E-02 

70 

GUNS 

NA 

9 

4499976E-02 

79 

.GUN  6 

NA 

9 

4499976E-02 

eo 

GWN7 

NA 

9 

4499976E-0P 

SI 

GWN3 

NA 

9 

4499976E-C2 

e? 

GUN9 

NA 

9 

4499976F-03 

63 

GW  10 

NA 

9. 

4499976E-02 

e4 

GUll 

NA 

9. 

449997  6r. -02 

85 

GUi2 

NA 

0 

OOOOOOOE-O 1 

86 

GEJ 

OBJ 

5 

00000002+01 

87 

STOREVN 

NA 

1 

43999996-00 

88 

STOREVNP 

NA 

-1 
c.  . 

605751CE-0) 

e9 

CR'JPN 

NA 

1 

6400000E-00 

90 

CRUPNP 

NA 

T 

40441238-01 

9] 

T07REVN 

NA 

T 

w 

cyryagcjiop  +  co 

9? 

TOTREVNP 

NA 

6 

393654 5E -01 

93 

APPMAX 

NA 

n 

w>  . 

5162797E-01 

OSECTION  P  —  VARIABLES 


VARIABLE 

INITIAL 

LOWER 

UPPER 

1 

NAME 

STATUS  VALUE 

LIMIT 

L  IMI  T 

1 

AREA 

5  OOOOOOCE+O 1 

0  OOOOOOOE-O 1 

NONE 

2 

STCR1 

3. OOOOOOOE+Ol 

5. 2990000E-0 1 

5 

50000002 +01 

3 

ST0R2 

3  OOOOCOOE  *-0 1 

5. 2990000E-01 

5 

SOOOOOOE+O 1 

4 

STQR3 

3  OOOOOOOE+Ol 

0.  2970000E-0 1 

5 

5COOOOOE  +  01 

5 

ST0R4 

3. OOOOOOOE+Ol 

5. 2990000E-01 

5. 

5000000E  +01 

6 

STOP  5 

3. OOOOOOOE+Ol 

5.  2990000E-0 1 

5 

5000000E+01 

7 

ST0R6 

3.  OOOOOOOE  *-01 

5.  2990000E-01 

5. 

50000002  *0 1 

26 


0 

ST0R7 

3  00G0C00E*01 

5 

299C000E-Q1 

5 

50000000 *01 

9 

STC'<3 

3  OOGOCCCO-ot 

5 

29900002-01 

5 

5COOOOOL *01 

10 

STCR9 

3.  0000000 L  ♦■0  l 

5 

89900008 -0 1 

5 

50000002 *01 

1  1 

STGR10 

3  OGOOCOOE  *  0 1 

5 

29900002-01 

5 

^0000002*01 

12 

STCRl  1 

3  0C0C0GCE*Ol 

5 

2990CG0E-0 1 

5 

50000002 *01 

13 

ST0R12 

3  00000000  *-01 

5 

2990000E-0 1 

5. 

50000002  *01 

14 

APP  l 

r>: 

o  ooooccco-ot 

0 

00000008-01 

0 

00000002-01 

10 

A  PP2 

FX 

0. OOOOCCOE-Ol 

0 

COOOOOOE-O 1 

0 

00000002-0 1 

10 

APP3 

FX 

0  00CGCGCC-01 

0 

OOOOOOOE-O 1 

0 

00000002-01 

17 

APP4 

1  89000002*0) 

0 

OOOOOOOE-Ol 

NONE 

10 

APP5 

1 . 89000000*01 

0 

00000002-01 

NONE 

19 

APP6 

1. 89000000*01 

0 

OOOOOOOE-Ol 

NONE 

20 

APP  7 

1.  Q9000C0C *0  1 

0 

OOOOOOOE-O 1 

NONE 

21 

APPG 

1. 89GOCCCC+01 

0 

OOOOOOOE-Ol 

NONE 

&.  e_ 

APP9 

1 . 89000000 *0 1 

0 

OOOOOOOE-O 1 

NONE 

23 

APP  10 

1 . 89000000*01 

0 

OOOOOOOE-Ol 

NONE 

24 

APP  1  1 

1 . 89000000*01 

0 

OOOOOOOE-Ol 

NONE 

25 

APP  12 

FX 

0  OOCOCOCO-Ol 

0 

OOOOOOCE-O 1 

0 

00000002-01 

26 

DRAIN  1 

LI. 

0  00000000-01 

0 

OOOOOOOE-Ol 

NONE 

27 

DRAIN? 

LL 

0. OOOOOCCO-Ol 

0. 

OOCOOOOE-O 1 

NONE 

23 

DRAIN3 

LL 

0  00000000*01 

0. 

OOOOOOOE-Ol 

NONE 

2? 

DRAIN4 

1 .  8900CC  01. +  0 1 

0 

OOOOOCOE-O i 

NONE 

30 

DRAINS 

1  39000000*01 

0 

00000002-0 i 

NONE 

31 

DRAINS 

1. 89000000*01 

0 

coccoooE-c i 

NONE 

* )  — 1 
«- 

DRAIN? 

i. e7ooocco*oi 

0 

0  OC 0000 E — G 1 

NONE 

33 

DRAINS 

1 . 09000000*01 

0 

OOOOOOOE-O 1 

NOME 

24 

DRAIN9 

1 . 39000000 *01 

0 

GOGCCOGE-O  1 

NONE 

33 

DRAIN  10 

1 . 3900CC CO *0 1 

0 

OOCOCOOE-C 1 

NONE 

36 

DRAIN  1  1 

1 . 39000000 *0 1 

0 

OOOOOOOE-Ol 

NONE 

37 

ORA IN12 

LI. 

0.  OOOOCCOC  *01 

0 

COOOOOOE-C 1 

NONE 

33 

ti  STCR1 

2  40000000-02 

0 

COOOOOCE-Ci 

NONE 

29 

N  5 TCP 3 

2.  4000CCCF.-02 

0 

OOOOOOOE-O 1 

NONE 

40 

N  STOP  3 

2.  4 000? P 00-02 

0 

OOOOOOOE-O i 

NONE 

41 

N  3 TOR 4 

2. 40000000-02 

0 

CCOOOCOE-O I 

NONE 

42 

N  3TCr«  5 

2. 40000000-02 

0. 

OOOOOOOE-O J 

NONE 

43 

N  3TCR6 

2.  40000000-02 

0 

OOOOOOOE-O 1 

NONE 

44 

N  STOP? 

2  40000000 -02 

0. 

OOC'OOOOE-C  1 

NONE 

40 

U  STCn3 

2.  40000C 00-02 

c 

OOOOOOOE-O) 

NONE 

46 

N  3 TCP  9 

2  400000.0-02 

0 

COOOOOOE-O J 

NONE 

47 

N  S ICR  10 

2  40000000-02 

0 

OOOOOOOE-O 1 

NONE 

43 

N  STCR  1  1 

2. 40000000-02 

0. 

OOOOOOOE-O 1 

NONE 

49 

N  STCR12 

2  4000C<  C  -02 

0 

OOOOCOOE-C 1 

NONE 

50 

CP OP  N! 

LL 

0.  00000 00!. -01 

0 

OOOOOOO E-0 1 

NONE 

51 

CROP  N? 

LL 

0  OOOOOOCC-OI 

0 

00000008-01 

NONE 

52 

CP OP  N3 

LL 

0.  00000008-01 

0 

OOOOOOOE-O  t 

NONE 

52 

CROP  N4 

5  0C0OCC0C-03 

0 

OOOOOOOE-Ol 

NONE 

S'l 

CROP  N5 

6  50CP0C;C-03 

0 

CCCOOOOE-O I 

NONE 

55 

CP  OP  NS 

5  OOOOOCC  C -03 

0 

OOOOOOOE-Ol 

NONE 

56 

CROP  N7 

6  5000G00C-03 

0 

00000002 -0 i 

NONE 

57 

CROP  N3 

6.  5000000L-03 

0. 

OOOOOOOE-Ol 

NOME 

53 

CROP  N9 

3  30000002-03 

0. 

OOOOOOOE-O) 

NONE 

59 

CROP  NIO 

LL 

o  oooo-cor.-oi 

0 

OOOOOOOE-Ol 

NONE 

60 

crop  N 1 1 

LL 

o  ooooooo t: -oi 

0 

OOOCOOGE-O 1 

NONE 

61 

CROP  N 1 2 

LI 

0  OOOOOGGL-OI 

0 

OOOOOOOE-Ol 

NONE 

27 


62 

N  PERC1 

LI. 

o.  owooooc::  -o  i 

0  0000000E-01 

NONE 

63 

N  PERCE 

LL 

0. 00000002-01 

0. 000000CE-G1 

NUNE 

64 

N  PEKC3 

LL 

0  OCOOCOCL-Ol 

0. GOOOOOOE-0 1 

NONE 

65 

N  P2RC4 

1. 890000CE-01 

0. OOOOOOOE-O 1 

NONE 

66 

N  PSRC5 

1.  89000002-01 

0. OOOOOOOE-O 1 

NONE 

67 

N  PERC6 

1 .  er/0G00CE-01 

0. OOOOCOOE-O 1 

NONE 

63 

N  PEKC7 

1 . Q900000E-01 

0.  OOOOOOOE-O  1 

NONE 

69 

N  PERCB 

1.  B9C-000CE-01 

0.  OOOOOOOE-O  1 

NONE 

70 

N  PERC9 

1. 89000002-01 

0. OOOOOOOE-O 1 

NONE 

71 

N  PERC10 

1. 8900000E-0 J 

0. OOOOOOOE-Ol 

NONE 

72 

N  PERC1 1 

1. 890000CE-01 

0. OOOOOOOE-O 1 

NONE 

73 

N  PERC12 

LL 

0, 000000 OL “01 

0. OOOOOOOE-O 1 

NONE 

74 

BYPASS  1 

FX 

0. OOOOOOCE-Ol 

0  OOOOOOOE-O 1 

0  00000002-01 

75 

BYPASS2 

FX 

0. OOOOOOOE -0 1 

0. OOOOOOOE-O 1 

0  00000002-01 

76 

BYPASS3 

FX 

0. OOOOOOCE-Ol 

0. OOOOOOOE-O 1 

0  00000002  01 

77 

BYPASS4 

FX 

0. 00000002-01 

0. OOOOOOOE-O 1 

0  OOOOOOOL-O 1 

73 

EYPASS5 

FX 

o. oooccces-oi 

0. OOOOOOOE-O 1 

o  00000002-01 

79 

BYPASS6 

FX 

0.  00000002-01 

0. OOOOOOOE-O 1 

0  00000002  01 

80 

BYPASS7 

FX 

0  00000002-01 

0. OOOOOOOE-O 1 

0  00000002-01 

81 

BYPAssa 

FX 

0.  00000002-01 

0. OOOOOOOE-O 1 

0. 00000002-01 

82 

BYPASS? 

FX 

0. 00000002-01 

0. OOOOOOOE-O 1 

0  00000002-01 

83 

BYPASS 10 

FX 

0.  00000002-01 

0. COOOOOOE-O 1 

0  00000002-01 

94 

BYPASS  1 1 

FX 

0  OCOOCCCE -0 1 

0. OOOOOOOE-O 1 

0  G0000C0L-01 

35 

BYPASS  12 

F  X 

0  00000002-01 

0. OCOOOOOE-Ol 

0. 00000002-01 

OTAL  FRACTICi 

,'AL  CHAN 

CE  IN  OBJECTIVE 

LESS  THAN  1  COOOCE-04  FDR 

FINAL 

RESULTS 

FOR 

3  CONSECUTIVE 

ITERATIONS 

0 

OSECTION  1  —  FUNCTIONS 


DISTANCE 


INITIAL 

FINAL 

FROM 

LAGRANGE 

NO 

NAME 

VALUE 

VALUE 

STX ~US 

NEAREST 

MULT 1PLIER 

BOUND 

1 

LAGOON  1 

0 

OOOOGE-O 1 

1. 

151002+01 

2 QUALITY  0 

0002  01 

o  oooooe-oi 

2 

LAG00N2 

0 

OOOOOE-O 1 

1 

151302  K>1 

EQUALITY  0 

OOOE  01 

0  OOOOOE-O 1 

3 

LAGQ0N3 

0 

000CCE-01 

1 

1  5 1  3C  t +0  ’ 

20UALITY  0 

CODE  0 1 

0  OOOOOE-O 1 

4 

LAGC3N4 

1 

64804E+0 1 

1. 

151 302  ►O  » 

EQUALITY-3 

41  It- 13 

0. OOOOOE-O 1 

5 

LAGOON 5 

1 

S9336E+01 

1. 

1 5130E +01 

EQUALITY-3 

41  IE-13 

0  OGOOOE-Ol 

6 

LAGC0N6 

1 

91 133E+01 

1 

15130E+0! 

EQUAL ITY-2 

160E  -09 

0. OOOOOE-O 1 

7 

LAG00N7 

l 

926V1E+01 

l 

151  30E  *-01 

EQUALITY-2 

356E-G9 

0  OOOOOE-O 1 

0 

LAG00N3 

1 

92609E+01 

1 

15130E  *-01 

EQUALITY-9 

840E-09 

0. OOOOOE-O 1 

9 

LAG0QN9 

l 

S9590E+0 1 

1 

151302+01 

EQUAL ITY-1 

Q90E  -08 

0  OOOOOE-O 1 

10 

LAGOON  10 

1 

3571 9E+0 1 

1 

151302+01 

EQUAL ITY-4 

4.22E  -03 

0  OOOOOE-O 1 

1 1 

LAGOONl 1 

1 

8366YE+01 

l 

1 5 130E +01 

EQUALITY-3 

649E-08 

0. OOOOOE-O 1 

12 

LAGOON  12 

0 

OGCOoE-O 1 

1 

151302+01 

EQUALITY  0 

coot  01 

0  OOOOOE-O 1 

13 

H2Q  BL1 

0 

OOOOOE-O 1 

0 

000002-0 1 

EQUALITY  0 

0001-  01 

0  OOOOOE-OI 

14 

H20  DL2 

0 

OOOOOE-O 1 

0 

000002-01 

EQUALITY  0 

coot  01 

0. OOOOOE-OI 

15 

H20  DL3 

0 

GOCOOE-O 1 

0 

OOCOOE-O 1 

EQUALITY  0 

OOOF  -01 

0. OOOOOE-O l 

16 

H20  DL4 

9 

03000E+00 

-6. 

670322-03 

EQUAL  IT Y-6 

670E-0S 

0  OOOOOE-O 1 

17 

H20  DL5 

-5 

5700 OE+OO 

-5 

252332-10 

EQUALITY-5 

252E-10 

0  OOOOOE-OI 

28 


13 

MCO 

f  L6  -6 

t  .-'OOCi  -OO 

-1 

lCcort-os 

EQUAL  I 7  Y- 1 

13/E-06 

0  OOOCOF-Ol 

1  ^ 

1 123 

2L7  -7 

9 "CO Of “00 

n 

3700  72  -."•fj 

EQUAL  I TY-3 

370E-03 

0  OOOOOE -0 1 

20 

1120 

DLo  -a 

37<;ocr  *cc* 

1 

43fct 1L-09 

EQUALITY  - 1 

466E-09 

0. OOOCOE-C 1 

21 

urn 

21.9  -5 

4 2000c ♦CO 

6 

8D304E  -09 

EQUAL 1TY-6 

E5GE- OR 

0, OOOOOE -0 1 

m 

c.  c_ 

H20 

DL 1 0-3 

6  70GGE  * 00 

-  1 

439 ! 2E -OS 

EQUALITY  -1 

439E--08 

0  OOOOOE-O 1 

m 

!  123 

CL1 1 -2 

•'2COOE  »00 

-3 

8O-67E-03 

EQUAL  I TY-3 

810E-09 

0  OOOOOE-O 1 

34 

!  120 

3L12  0 

OOOOOE-O 1 

0 

OOOOOE-O 1 

EQUALITY  0 

OOOE  01 

0. OOOOOE-O 1 

25 

dra 

CP  1  4 

jf'i'TF-O  1 

1 

1  167C.L  t-02 

FREE  1 

1  i  /E  ►or 

L 

36 

OH  A 

CP2  4 

3?<,G«"t>01 

1 

1  167c-2  ♦02 

FREE  I 

1 17E 702. 

L 

27 

DR  A 

CP3  4 

3 OOOGE+O 1 

1 

1  1  67-12  ♦03 

FREE  1 

1  17E<-02 

L 

23 

DHA 

CP4  2 

4iOO0E«01 

6 

2644  71.+0) 

FREE  6. 

264F+01 

L 

29 

DR  A 

CP5  2 

4  6000E+0 1 

1 

033C9C-OD 

FREE  ! 

033F+02. 

L 

30 

DR  A 

CP6  2 

460CCE+0 1 

1 

022 1 iE ♦Or 

FREE  1 

C22E  *02 

L 

31 

DRA 

CP7  2 

46000£«  0 1 

1 

023B2C  *-02 

FREE  1 

024E702 

L 

32 

DR  A 

CPS  2 

460C0E*  0 1 

a 

B7446E+01 

FREE  e 

874E+G1 

L 

33 

DRA 

CP9  2 

46000E+01 

9 

7660  7t.  ♦Ol 

FREE  9 

766E+01 

L 

34 

DRA 

CP10  2 

46000E+0 1 

1 

06542C«-02 

FREE  ! 

C6DE+02. 

L 

35 

DRA 

CPI  1  2 

S?COCE»CO 

3. 

421 CtE  ♦O 1 

FREE  3 

42  IF.  +01 

L 

36 

DRA 

CP  1 2  4 

35000E»0 1 

1 

1  1O76E+02 

FREE  1 

1  1  7E  t-02 

L 

37 

N  CAL  1  0. 

OOOCOE-O 1 

0 

000007.  -0 1 

EQUALITY  0 

OOOE  01 

0  OOOOOE-O 1 

33 

N  BAL2  0 

00000E-01 

0 

GOOGCE-O 1 

EQUALITY  0 

OOOE  01 

0. OOOOOE-O 1 

39 

N  BAL3  0 

OOOCOE-O 1  0  000002 -0 ) 

FINAL  RESULTS 

EQUALITY  0 

OOOE  -01 

0  OOOOOE-O 1 

SECTION  1  -- 

NO.  NAME 

FUNCTIONS 

INITIAL 

VALUE 

FINAL 

VALUE 

STATUS 

DISTANCE 

DROM 

NEAREST 

BOUND 

LAC-RANGE 
MULT  I  PL  I  EH 

10 

N 

CAL4 

-l 

46C00E-02-6 

302000-05 

EQUALITY -6 

332E-05 

0  OOOOOE-O 1 

4  1 

N 

CAL  5 

o 

04000E -02- 1 

831=90-05 

EQUAL  I  TV- 1 

S32E--05 

0.  OOOOOE-O 1 

42 

N 

3AL6 

-1. 

46000E-02-6 

408710-05 

EQUAL 1 TY -6 

409E-05 

0. OOOOOE-O l 

43 

N 

BAL7 

6 

04000E-02-2 

720100-05 

EQUALITY -2 

720E-05 

0. OOOOOE-O 1 

44 

N 

3AL8 

6 

0400QE-09-3 

402600  -05 

EQUAL  1 TY-3 

4030-05 

0  OOOOOE-O 1 

4  5 

N 

BAL9 

-9 

96000E-C2  1 

67  5 1 5  0  -05 

hOUALITY-1 

67DE-05 

0. OOOOOE-O 1 

4w» 

N 

BAL10 

_ 

CL  . 

64600E-01  1 

45S31E-13 

EQUALITY- 1 

4  58E-13 

0  OOOOOE-O 1 

47 

N 

BALI  1 

_ 

e. . 

64600E-0 1 -3. 

52306E-06 

EQUAL I TY-3 

523E-06 

0. OOOOOE-O 1 

40 

N 

BAL12 

0 

OCOCOE-O 1  0 

000000-01 

EQUALITY  0 

OOOE-  01 

0  OOOOOE-O 1 

49 

N 

CROP  l 

0 

OCOCCE-C-l  *'• 

0005-2-01 

EQUALITY  0 

' OOF  Cl 

0, OOOOOE-O 1 

50 

N 

CROPS 

0 

OOOOCE-Ol  0 

000002-01 

EQUALITY  0 

.000  01 

0. OOOOOE-O 1 

51 

N 

CRQP3 

0 

OOOOCE-Ol  0 

OCOCOO-O t 

EQUALITY  0 

0000  01 

0  OOOOOE -01 

52 

N 

CROP  4 

-1 

03 1 37E-04  5 

01093C-05 

EQUALITY  -5 

01  10-05 

0  OOOOOE-O 1 

53 

N 

CR0P5 

6 

808olE-04-l 

276950-07 

EQUALITY-1 

277E- 07 

0.  OOOOOE-O 1 

54 

N 

CR0P6 

-1 

03137E-04  6 

180130-05 

EQUAL  I TY-6 

1  SUE-05 

0  OOOOOE-O 1 

55 

N 

CP0P7 

6 

S036 1 E-04 - 1 

64 62 rE -07 

EQUAL 1TY-1 

6-46L-  07 

0  OOOOOE-O 1 

56 

N 

CP0P8 

6 

8036 1 £-04- 1 

69s34t— 7 

EQUAL  I T Y-  1 

6966—07 

0.  OOOOOE-O  t 

57 

N 

CR0P9 

-7 

1  7948E-04  '5 

1 E372E-03 

CGUALITY-5 

1S9E-08 

0.  OOOOOE-C 1 

53 

N 

CROP  1 0 

0 

OOOOOE-O I  1. 

l 3566E  -  l  5 

EQUAL  I T Y- l 

136E-15 

0.  OOOOOE-O 1 

59 

N 

CROP  1 1 

0 

OOOOCE-Ol  5 

00 1 1 OE - 1 5 

EQUALITY-5 

00 1 E  - 1 5 

0  OOOOOE-O 1 

60 

N 

CROP  12 

0 

OOOOOE -01  0 

000000-01 

EQUALITY  0 

CGOE  01 

0  OOOOOE-O 1 

61 

LAG  N 1 

7 

2000CE-02  4 

605000-01 

EQUALITY-9 

343E-1 1 

0  OOOOOE-O l 

29 


<  n 

Wv- 

LAG  N2 

7  2C00GE-02 

4  6060CE-OI 

OGUAL I TY 

-4 

6273  1 1 

63 

LAG  N3 

7  2G000E-G2 

4. 60SOOC-OJ 

LGUALITY 

-6 

774L-1  1 

64 

LAG  N4 

5  97toor.-oi 

a  ioic -oi 

i'ou^lj  ry 

-7 

4 SOL- 1 1 

65 

LAG  N5 

5.  976C0E-01 

4  605O0E-01 

EQUAL  I TY 

-4 

2C6E-1  1 

66 

LAG  N6 

5  97600E-01 

4, 60500L-01 

i:  DUALITY 

CL  . 

1G0L  -08 

67 

LAG  N7 

5  R7600E-01 

4  60500C-0 1 

rOUALITY 

-5 

077E  08 

60 

LAG  N3 

5  97600E-01 

4  60 4 ^7 L -01 

EQUAL  1 TY 

-3 

29 1  E-Os 

69 

LAG  N9 

5  R^QOE-Ol 

4.  60400:.  -01 

EQUALITY 

-9 

RP0E-G5 

70 

LAG  N 10 

5  R7600E-01 

4  605-SC-Q 1 

EQUALITY 

-5 

94  SE  05 

71 

LAG  Nil 

5  R7600E-01 

4  6030011-01 

EQUALITY 

-1. 

593C-07 

72 

LAG  N12 

7  200G0E-02 

4. 6050GC -0 1 

EQUALITY 

-1 

1 25F  -10 

73 

GW  GUAL 

-1  421  ORE- 14- 

-2  2S3o7E-09 

LCLLRDND 

_n 

2S4E-0-?  L 

74 

GWN1 

0  OOOOOE-O 1 

0. OOOOCE-O 1 

IGNORED 

i 

OOOF  EJ3 

75 

GWN2 

0  OOOOOE-O 1 

0  OOOOOE-O 1 

IGNORED 

i 

OCOE  --30 

76 

GWN3 

0  OOOCOE-O 1 

0.  OOCOOE-Ol 

IGNORED 

i 

COOL*  • 30 

77 

GWN4 

9, 4500CE-02 

6  4ii  IDE-Gl 

IGNORED 

i 

OOOL"  ■■■20 

70 

GWNO 

9  450C0E-02 

9. 37A3_G-02 

I  GNCRED 

i 

OOOF -3A 

79 

GWN6 

9. 45000E-02 

6. 2QC36C-02 

IGNORED 

i 

OCOE --3  0 

1 

FINAL  RESULTS 

0 

OSECT I  ON  1  — 

FUNCTIONS 

D I  STANCE 

IN IT IAL 

FINAL 

IRON 

NO 

NAME 

VAL  UE 

VALUE 

STATUS 

NEARED r 
BOUND 

80 

GWN7 

9. 450C0E-02 

4. 940070-02 

IGNORED 

i 

OOOF  ■  :•  - 

81 

GWN3 

9.  45000E-02 

7  519400-02 

IGNORED 

i 

OOOF : 3 

82 

GWN9 

9  45000E-02 

-2.  230900-03 

IGNORED 

i 

OCOE f .  ■ 

83 

GW  10 

9.  45O00E-02 

7  70 1 31 C-02 

IGNORED 

1 

r$c  ( »  -  t  3 

84 

GW1 1 

9.  45000E-02- 

-3. 030360-01 

IGNORED 

l. 

00 Oe  •  Iv 

35 

GW  1 2 

0. OOOOOE-O 1 

0  OOCOOE-Ol 

IGNORED 

i 

CODE s 3 7 

86 

OBJ 

5  OOOOGE+O 1 

1  203630  *-02 

OCR 

37 

STOREVN 

1  44000E+00 

1 .  73~  2^2 +G0 

IGNORED 

i 

07  K‘[- J-  27' 

83 

STORE VNP 

2  60575E-01 

3  1 47320-01 

I  i> 

i 

OC-LE-  " 

39 

CRUFN 

I.  oJOOOE-'CO 

2  385770  m'-o 

ICNOREO 

1 

. 

=50 

CHUFNP 

3-  404  4  IE-01 

4  v«-r.-r.L-nj 

IGNORED 

i 

v7-  i-  • 

=U 

TCTREVfJ 

3  OSOCCE^OO 

.1  124  76L  »00 

IvN  jRE  iv 

i 

'Or  '■  . 

92 

TOTRCVNP 

6  3^3655-01 

3  D6c.*i  1 

IGNORED 

i 

OCOE ■ 30 

°3 

APPMAX 

3  5162LE-01 

1  36.3 1 SC-0 1 

IGNORED 

i 

OOC'L  <  30 

1 

FINAL  RESULTS 

0 


OSECT ION  2  — 

-  VARIABLES 

NO  NAME 

INITIAL 

VALUE 

FINAL 

VALUr 

STATUS 

D  I  ST  ANTE 
!  NOE' 

near;.  ■ 
BOUND 

1  AREA 

5.  OOOOOE+0 1 

1. 20363E+02 

SUPBASIC 

1  2S4E*03  i 

0  00000E-01 
0. 0000GE-01 
0  OOOOOE-O 1 
0  OOOOOE-O 1 
0. OOOOOE-O 1 
0  OOOOOE-O 1 
0  OOOOOE-O 1 
O. OOOOOE-O 1 
0  OOOOOE-O 1 
0  OOOOOE-O 1 
0  OOOOOE-O 1 
0  OOOOOE-O 1 


L  At'K 

mui.  i  i flier 


reduced 

<, : < Ai <  j  snt 


1  OOuOUL*GO 


30 


n 

GTCR  1 

-> 

c:o"<-:e+oi 

3  1970011+01 

BASIC 

303r*-01  u 

3 

21C22 

n 

CCOOCE  *-C  1 

4  34B70E+01 

BASIC 

i 

151E*-01  U 

o 

ST  09  3 

n 

CC0002+0J 

0  300002+01 

NCN3ASIC 

UPP2RHND 

0  OOOOOE-Ol 

5 

STCR4 

-} 

OC  uOGE  -*-0 1 

5  50G00E+01 

NOun as i c 

UPPLKUND 

0.  OOOOOE-Ol 

6 

STCRS 

"J 

000002+01 

5  5CC002+0 1 

NGNBASIC 

UPP2R/IND 

0.  OOOOOE-Ol 

7 

STOPS 

-5 

OOOOCE+O 1 

0  732922+01 

DASIC 

7 

1 5 1 E  +  00  U 

3 

STGR7 

3 

OOOOCE+O 1 

3  714632+01 

BASIC 

1 

7S5F+01  U 

7 

STOPO 

OOOOCE+O 1 

1  023712+01 

BASIC 

1 

031E  +  01  .  L 

10 

ST0R9 

n 

w* 

OOOOOE+Ol 

3  320012+00 

BASIC 

" 1 
G. 

790F+00  L 

1  1 

ST3R 10 

3 

OCOOOE+O 1 

1  54315E+01 

BASIC 

1 

49CE+01 : L 

12 

STOP  1 1 

3 

OOOOCE+O 1 

8  94  EOGE+OO 

BASIC 

8 

4 i 82+00  L 

13 

ST0R12 

n 

CCOOCE +01 

2  04  c?  1  Ot+01 

BASIC 

1 

993H  +  01  :  L 

10 

APP  1 

0 

OCOOOE-O 1 

0  OOOOOE-Ol 

NONBASIC 

FIXED 

0. OOOOOE-Ol 

13 

APP  2 

0 

OOOOOE-Ol 

0  000002-01 

NONBASIC 

FIXED 

0. OOOOOE-Ol 

16 

APP3 

0 

OOOOOE-O 1 

0. OOOOCC-Ol 

NOMBrtSIC 

FIXED 

0. OOOOOE-Ol 

17 

APP4 

1, 

8P00CE+01 

1.  5949CL  +0 1 

BASIC 

1 

595E+01 : L 

18 

APP  5 

1 

Q^OOOE+O 1 

1  15431E+01 

BASIC 

1 

154E+01  2 

19 

APP  6 

1 

89 OOOE+O 1 

1.  82327E  t-0 1 

BASIC 

1 

82GE+G1  L 

20 

APP  7 

1 

5P000E+01 

2  13294E+01 

SUPSASIC 

tC 

.  1532 i 01  .  L 

0. OOOOOE-Ol 

21 

APP  3 

1. 

S~O0CE+O 1 

3.  71603E+01 

SUPBASIC 

3 

716E  +  0I  :  L 

0.  OOOOOE-O  i 

H  c. 

APP  9 

1, 

89000E+01 

1 . 9 10S2E+01 

SUP BASIC 

1 

91  IE  i-O!  L 

0. OOOCOE+Ol 

23 

AFP  10 

1 

09 OOOE+O 1 

-4  o 4  4  0 e-C  - 1  4 

NuNpASIC 

LCUERilMD 

0  OOOOOE-Ol 

20 

APr  1  1 

1 

ecoov-E+o  i 

1  8970-2+01 

SUPBASIC 

I 

8902+01  L 

0  OOOOOE-O 1 

23 

APP  '  2 

0 

OOOOOE-Ol 

0.  0-' >0002-01 

NONBASIC 

FIXED 

0. OOOOOE-Ol 

26 

DR A  INI 

0 

OvOOCG-O 1 

0. 00CGGL-01 

BASIC 

LGU2RBND 

t; 

c_  / 

DR A I M2 

G 

OOOOOE-O 1 

0. 000002-01 

BASIC 

LOWERLND 

23 

DRAIN3 

0. 

OOOCOE-O 1 

0. 000002-01 

BASIC 

LONER HMD 

29 

DRAINO 

1 

39000E+01 

4. 903152+01 

DASIC 

4 

90Gr  !-0 1  .  L 

20 

DRAINS 

1. 

S^OOGc+O 1 

8.  336S3L+00 

BASIC 

8 

337E+00  L 

31 

DRAIN6 

1 

89000E+01 

9.  46047C+00 

BASIC 

9 

460E+00. L 

32 

DRAIN7 

1 

3 '’OOOE  +  O  1 

9  294202+00 

BAS  1 0 

9 

294E+C0: L 

n 

DRA  r  M3 

1 

s^cooe+oi 

2  29316E+01 

BASIC 

2 

2932+01  L 

30 

DRAIN9 

1 

S9000E+0 1 

1  40 1 55L+01 

BASIC 

1 

4022 <  01  :  L 

23 

DRAIN10 

1 

S^OCOC+Ol 

5  134542+00 

BASIC 

5 

1352+00.  L 

36 

DRA INI  1 

1 

8  ‘'J  OOOE+O  1 

2.  16275E+01 

DASIC 

2 

1 63E+0 1 : L 

3  V 

DRAIN12 

0 

0CCC~E-0 1 

0  000002-01 

BASIC 

LCwrrusND 

-  n 

N  5 TOR  1 

■1CC0CE-02 

3. 32774E-02 

BASIC 

3 

32GE-02  L 

29 

H  5  T  C  R  2 

c. 

•■•000CE-O2 

3  222232-02 

BASIC 

n 

2222-02  L 

FINAL  RESULTS 

SECTION  2  — 

VARIABLES 

DISTANCE 

INITIAL 

F I NAL 

FROM 

REDUCED 

NO 

NAME 

VALUE 

VALUE 

STATUS 

NEAREST 

GRADIENT 

BOUND 

40 

N 

ST0R3 

2 

40000E-02 

3 

1 024°E-02 

BASIC 

3 

102E-02 

L 

41 

N 

3  TOR  4 

*- 

“000CE-02 

2 

553962-02 

BASIC 

n 

554E-02. 

L 

42 

N 

STCR5 

T 

e. 

4COOOC-C2 

2. 

3300SE-02 

BASIC 

o 

. 

380E-02 : 

L 

43 

N 

STCR6 

2 

40C0CE-02 

c 

;’9995i:-02 

BASIC 

c.  . 

300E-02: 

L 

44 

N 

STCR7 

e. 

4C00CE-02 

n 
e.  - 

331 40E-02 

BASIC 

c~  . 

331E-02 

L 

31 


45 

N  STOPS 

n 

c_ 

* C^OCC-rp 

2  8002! E- 02 

BASIC 

2  902E  02 

(. 

46 

N  STOP 9 

n 
e.  . 

4CGCCE-02 

3.  9327 IE -02 

BASIC 

3.  933E-02. 

L 

47 

N  STCR10 

n 

C- 

4C000E-02 

3.  412622-02 

BASIC 

2  4  ME -02 

L 

43 

N  STOP  1  1 

A  . 

4G0C0E-02 

3  12143C-G2 

BASIC 

3  121E-02 

L 

47 

N  S  TOR  1 3 

-> 

c. 

4000CE-G2 

3.  373512-02 

BASIC 

3  379E-02 

L 

50 

CRC?  N! 

0 

OOOOOE-Ol 

0.  000002-01 

BASIC 

LOWERBND 

51 

CROP  NO 

0 

OCCCCC-OJ 

0  00C0CL-01 

BASIC 

LCW2RDND 

50 

CROP  N3 

0 

00000E-01 

0  0000CL-01 

BASIC 

LCW2RDND 

53 

CROP  N4 

5 

C000CE-03 

2  31741L-03 

BAS  1 C 

2  317E-C3 

L 

54 

CROP  N5 

6 

5000GE-03 

l  975992 -03 

BASIC 

1  976E-03 

L 

55 

CROP  N6 

5 

■OGOOGE-03 

2. 72C1GC-03 

BASIC 

2  720E-03 

L 

56 

CROP  N7 

6 

50CC0E-03 

3. 17914E-03 

BASIC 

3  179E-03 

L 

57 

CROP  NS 

6 

5G000E-03 

4. 96294E-03 

BASIC 

4  963E-G3 

1. 

58 

CROP  N9 

3 

3G0CGE-G3 

2  9273cE-03 

BASIC 

2  9272-03 

L 

59 

CROP  N 1 0 

0 

OOOOOE-C 1 

1  13566E-15 

BASIC 

LOWERBND 

60 

CROP  Nil 

0 

OGGCOE-O 1 

5  001  ICE -  I 5 

BASIC 

LOWERBND 

61 

CROP  N12 

0 

OOOGOE-O 1 

0. OOOOOE-Ol 

BASIC 

LOWERBND 

62 

N  PERC1 

0 

OOOCL  Zi-0 1 

0  000 OC  E - 0 1 

NGN BAS ! C 

0  000C02-0 1 

63 

N  PERC2 

0 

CCCOOE-O 1 

0.  OCOCCE-O 1 

NCN3AS1C 

LOWER BN D 

0  OCOCK  2-01 

64 

N  P ERC3 

0 

GGOOOE-O 1 

0.  OOOOCE  -01 

H0N3ASIC 

LOWERBND 

0. OOOOCE *01 

65 

N  PERC4 

1 

39C00E-0 1 

O  qor.-cr-.pp 

BASIC 

0  936E-02 

L 

N  PEKC5 

1 

3" GOGS -0 1 

2.  U07?;;-02 

BASIC 

3  11)2-03 

L 

6V 

N  PERC6 

l 

39GG0E-0 1 

7.  9 1 0332-02 

BAS  I  C 

7  9 ioE-02 

1. 

63 

N  PERC7 

L 

9900CE-0I 

9  001232-02 

BASIC 

9  0012-02 

L 

69 

N  FERC8 

1 

o’-CC  C  F.  -  0 1 

2.  637 -6-2-01 

BAS  I C 

2  653E-01 

L 

70 

N  PER 79 

L 

390COE-0) 

2  1 3"-  142-01 

BASIC 

2  1251-1-01 

L 

7) 

N  PERCiO 

I 

8900CE-01 

-1.  576202- 1  5 

NON 3 AS  I C 

LOWERBND 

0  OOOOOE-C 1 

-r—i 

r  r- 

N  PERC 1 1 

1 

5^ 00 OF -0  1 

6  2  7  4  ^  .--i:  “0  J 

BAS  1 C 

6  274E-01 

L 

73 

N  PERC 12 

0 

OCOCCE-O 1 

0. 000002-01 

NGNDASIC 

LOWER;  N“ 

0  OOOOOE-Ol 

74 

BYPASS  1 

0 

OvOO<.’i--v  1 

0  Ov^  -■  -•-»  1 

NONDAS10 

:T  ?  v  • :  . , 

0  OOOOOE-C' i 

75 

BVPASS2 

0 

OOOCOE-O 1 

0  OOOOCE -01 

N0NDA5IC 

FIXEB 

0  OOOOOE-Ol 

7cj 

3  Yp A533 

0 

OOGCCE-O 1 

o.  ooc-.  :t:  -01 

MONO AS) C 

FIX  f-  i  '• 

0  OOOOOE-C l 

7  7 

BYPASS 4 

0 

OOCCOE-C 1 

0.  000‘ <02  -0  1 

NON BAS  1C 

FUEL 

0.  COOOOE  -0 i 

73 

BYPf  .SS5 

G 

OOOOOE-O 1 

0. OOOOOE-Ol 

NO;  JB  AS  1C 

r  i  a  u  i  > 

0.  OOOOOE-O  1 

79 

3YPASS6 

0 

OOOCOE-O 1 

0  000002-01 

NGNDASIC 

FIXED 

0. OOOOOE-G l 

FINAL  RESULTS 


OSECTION 


VARIABLES 


DISTANCE 


INITIAL 

FINAl 

FROM 

REDUCED 

JO 

NAME 

VALUE 

value 

STATUS 

NEAREST 

GRADIENT 

BOUND 

30 

3YPASS7 

0 

OOOOOE-C 1 

0  COOOOE -01 

NGNDASIC 

FIXED 

0 

000002-0 1 

81 

3YPASS8 

0 

OOOOCE -0 1 

0  003002-01 

NGNDASIC 

FIXED 

0 

OOOOOE-O 1 

82 

BYPASS9 

0 

OOOOOE-Ol 

0. 000002-01 

NONE ASIC 

FIXED 

0 

OOOOOE-O 1 

S3 

BYPASS10 

0 

COOOOE-O 1 

0  OOOOOE-Ol 

NGNDASIC 

F  VGA 

0 

OOOOOE-Ol 

34 

BYPASS l 1 

0 

OOCOOE-O 1 

0  OOOOCE -Ql 

NGNDASIC 

F ;  xi- 1) 

0 

OOOOOE-O 1 

es 

BYPASS  12 

0 

GGCOGE— 0 1 

0.  OOOOOE-Ol 

NuNBASIC 

F  i  xED 

0 

OOOOGE-O l 

1  RUN  STATISTICS 
0 


ONUMBER  OF  ONE  -DIMENSIONAL  SEARCHES  =  71 

ONEWION  CALLS  »  304  NEWTON  ITERATIONS  =  267  AVERAGE  = 

0FUNCT10N  CALLS  =  569  GRADIENT  CALLS  =  70 

ACTUAL  FUNCTION  CALLS  (INC  FOR  GRADIENT)  =  6689 

ONUMBER  OF  TIMES  BASIC  VARIABLE  VIOLATED  A  BOUND  =  44 

ONUMBER  OF  TIMES  NEWTON  FAILED  TO  CONVERGE  «-  17 

OTINGS  STCPSIZE  CUT  BACK  DUE  TO  NEWTON  FAILURE  =  14 


32 


APPENDIX  B:  LTMOD  program  listing 


C  LIMQD  PROGRAM 
C 

C  MAXMEM  SETS  Thh.  i'.MC'JJ  T  OF  MEMORY  USED 
PARAMETER  I-1AXMEM  =  20000 
* 1 NSERT  SY5C DM:  AH-  GYP 
REAL»0  2  <  M,V*  ME  M  ) 

COMMON  /FUPCE/Z 
INTEGER  *4  NC0RE 

INTEGER  *2  INFILL  3).  GUTFIL<6>,  I  NAME  (  1  6  » ,  OUTNAM  (  X  e.  ) 

I  NTECER  -*2  CGN6G!5  > .  CGNAME  !  16  ) 

INTEGER  *2  DEFILE  16).  DEE NAFF  <  16) 

LCCICAL  LOG 

DATA  NCORE/UAXMEM/ 

DATA  INF ILE,  ) CM  INPUT  FiLC/ 

DATA  CU1FIL  '  1  liiOL'TPUT  FILE/ 

DATA  CCMSC,' lOt-.r^/.M  FILE/ 

DATA  DEFILE/ 11  ml EFINE  FILE  / 

WRITE! 1  , 100) 

100  FORMAT (  'LAND  TREATMENT  OPTIMIZATION  MODEL. '» 

10  L0G=0"VP-*A<  INFIi_£  INTO! 10  « ,  ALREAD,  I  NAME,  I'NTS!  32)  .  INTSi  1>. 

*  AiVOLD.  0.0) 

IF  {.NOT  LOO)  C-0  TO  10 

15  LOG-OP VP*A  ( CCMSO  •  i NTS ( 1 0 > , a IREAD- CONAME.  INTE'02)-  '  N  TEC*), 

*  ALVOLD,  0,  O) 

IF ( .  NOT  LOO)  GO  TO  15 

16  LOO-OP  VP -vA  <  DEFILE.  INTO  (  1 0  ) ,  AtRCAD.  DEFNAME,  INTO  <  32  •  •  I N  1 0  <  0  1 

*  ALVOLD, 0. 0) 

I F  (  MOT.  LOO).  CO  TO  16 

20  LOG-OP VP  *A!0U7FIL,  INT5( 11 ) ,  ASWRIT . OUT  NAN,  INTO ( 32  I  NT 0<  ?) . 

*  ALOVAP, 0, 0 ) 

IF<  NOT.  LOG)  GO  TO  20 
CALL  GRG(Z. NCORE) 

LOG  =  CLOSLA!  INTSd)) 

LOG  -  CLOSiA ( into; 2) ) 

LOC  =  CLOSiA! INTO- 3) ) 

LOG  =  CLOG-LA!  INTS<4)  ) 

CALL  EXIT 
END 
C 

SUBROUTINE  GCCMF  !G.  X) 

IMPLICIT  REAL  •  A- 7  ) 

C  OMMON / I N 1 T  DH / I N I T 

REAL  H12>,P(i2),E(12),W<;2).  MDR  (  1?) ,  NI  (  12  )  . 

*  NK<  12  )  ,  DEN!  12)  .  CP!  12),  MNU  i  12).  FNL  !  12  ) ,  NCOM  < 1 2  > • EFF ( 1 2  > 
SAVE  I. P. £, W, MDR, NI.  NK,  DEN,  CP,  MNU.  FNL,  NCOM.  EFF , N~  ' 

*  OA,  OB. DP 

IF! IN I T.  EG  0)  GO  TO  30 

READ! 7.  10) I. P. E. W, MDR.  NI  .  Ni\.  DEN,  CP,  MNU,  FNL,  NCOM, EFF, NCO, 

*  OA,  03. DP 

10  FORMAT ! 26  <  6F 1 0  6/1.4F10  6) 

WR  l.TE  (  6,  12) 

12  FORMAT!  'PARAMETERS ' ) 

WRITE! 6.  1  1  )  I . P. E, W, MDR, NI,  NK,  DEN, CP, MNU, FNL,  NCOM,  Erl  , NCO. 

*  OA, 03, DP 
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f  1. 1  MOD  PROGRAM 


11  FORMAT! 26  (  IX,  tno  6/>.  IX,  4F  10.  6) 

30  CONTINUE 

DIMENSION  0(93). X (05) 

C 

C  STORAGE  WATEf?  BALANCE 
C 

G(  1  )=X  <  14>+X(74)-t-X<2)-Xi  13)  +  <E<  1)-P<  1  )  >  e-DriAXl  (  X  (I 

*  X ( 4  ) .  X  <  5  > .  X ( 6) .  X  <  7  ) .  X(S).  X  (9) ,  X(10),  X(ll  ),  '<  <  12)  . 
G  (2  )*X  (  1 5)  «•*  <  75>  *■<  (  3)  -X  <2  >  +  <E<  2  i-P  (2)  )*DMAX  1  (  X  (2  . 

*  X<4)»  X  (  5  > .  X  ( 6  >  •  X(7).  X  ( 8  > .  X(9).  X  (  1 0 ' .  X  (  1 1  )  <  X  <  1 2 )  ■ 
G(3)=X  <  16  >>X  (  76)  t-X  <4)-X<3>  +  <E<3>-P  (  3>  >  *DMAX1  (  X  (2  ! 

*  X ( 4  ) .  X ( 5  )  .  X  (  6  > .  X  <  7  ) .  X  <  3  )  .  X<9),  X  (  1 0  > .  X  (  1  I  ),  X( 12), 
C(4)=X( 17 ) +X ( 77 )  +  X  <  5) -X ( 4 )  +  <  E  <  4 ) -P ( 4 ) ) *DMAX 1 ( X  <  2 1 

*  X(4),  X  (  5  )  .  X  <  6  >  •  X  <  7  ) ,  X  <  3  > ,  X  (  9 ) ,  X  <  1 0 ) .  X(ll).  X  <  1  2 ) , 
G(5)*X  <  13)  *-X<73)+X(6)-X<5)  +  (E<  5)  -P  <  5)  )  *DMAX1  (  X  (2: 

*  X  <  4  ) ,  X  (  5  ■  .  <  ( 6 ) ,  X  ( 7 ) ,  X  ( e ) .  X  ( 9 ) .  X  (  1 0 ) .  X  ( 1  1  ) .  X  (  1 2  > , 
G  <  6 )  =  X  <  1 9  )  *■  x  ( 79 )  ■►  x  (  7  >  - X  ( o )  -1-  ( E  <  v>  >  --P  <  6  >  )  -  DM AX  1  (  X  (2  : 

*  X  (  4  ) ,  X  (  5  i ,  x  ( 6 ) ,  « i  7  ) ,  X  ( 3  > .  X  <  •>  > ,  X  <  1 0  ) »  X  (  1  1  ) ,  X  <  12). 
G(7)=X (20>  rS <8C)+X (3)-X (7 '  ► (£(7)  -P<7> ) *DMAXi < X (2 : 

*X(4),  X  <  j  )  .  X  {  6  ) .  X  (  7  ) .  X  (  3 )  •  X  !  9  ) ,  X  <  1  0  > ,  Kill  >,  X  i  12.').  ) 
G<8)  =  x  (21  )*X  (31  >+X  (O)-M9w(E(0)-P  (B>  >  *DMAX  1  <  X  <2  ! 

*  X  (  4  ) .  X  i  5  >  .  X  <  6  ) ,  X  i  7  ) ,  X  (  3  )  ■  *  (  9  > ,  X  <  1 0  )  <  X  (  1 1  )  .  X  (  1  2. )  « 
G  ( 9 )  =  X  <  22  )  ♦  X  <  32  '  + «  <  1 0 )  - X  <  9  •,  +  <  e  (  9  >  -P  i  9  )  )  ►DMA  X  1  (  X  i  Z 

*  X  (  4  ) ,  X  (  5  '  •  X  (  6  )  .  X  (  7  i ,  X  ( P  )  ,  x  (  Q  ) ,  \  {  i  o  > ,  X  <  1  1  .  X  f  1  2 )  . 

G  (  1 0  >  =  X  (  22  ■  -  *  •  33  ►  X  i  1  1  )  - X  >  10  )  ►  ( E  (  10 )  ~P  (  1 0 )  )  ►DMA-  1 

*  X  (  4  ) ,  X  i  5  >  •  X  (  3  ■ X  >’  7  :• ,  X  i  8  )  •  «  ( 9  ) ,  X  <  1 0  ) .  X  <  1  1  > ,  X  (1 2  ) 
G< 1 1 >  =  X  <24 '♦X . 34 itx < 12  >-< •  11  >  +  <E( 11 >-P( 1 1 > )*9MAV 1 

*  X  ( 4  > ,  X  <  5  > .  X  ( 6  > ,  X  v  7  > ,  X  <  3  > .  X  <  9 ) ,  X  <  1 0 ) .  X  <  1 1  ) .  X  U  2  >  • 
C  <  12  >  =X  (  25  1  ♦X  <  8S )  ■>■  X  <  1 3  '  -  X  <  1  2 )  +  <  E  (  12) -P(  12)  )*DMAX1 

*  X  (  4  ) ,  X  <  5  ) .  X  <  b  > ,  X  (  7 ) ,  X  ( 3 ) .  x  ( 9 ) ,  x  (  1 0 ) .  X  (  1 1  )  ■  <  <  1  2  )  . 
C 

C  WATER  BALANCE  AT  IRRIGATION  SITE 
C 

G( 13)=EFF( 1 ) v<  X  < 14 ) rx  <74 ) )  ♦  (P ( 1  >-W< 1 ) >  *X ( 1  ) -x (26 , 
C  (1 4  )  =EFF  <  2  >  *<  X  (  1  5  >  *X  <  75  >  >  +  <P  <  I? )  -W  i  ?  >  >  *X  <  1  >  -X  <  2?  ' 
G  v  15 )  =EFF  <  3  )  * \  X(  i6>  *-X<76)  i  +  (PCD  -W<3>  )*X(  1  )  -X  (20. 
G<  16  )  =EFF  <  4  ■  *  (  X  (  1  7  )  *X  (  77  >  1  +  <.P  (  4  >  -W(  4  )  >  »X  <  1  )  -X  (  2«  • 
G  >  1 7  )  =FFF  (  5  )  v  (  X  (  1  ?, '  *-X  (  73  )  )  ♦  <  P  ( t, )  -W  (  5 )  >  »X  v  1  )  -  X  (  30  ; 
G  <  10 )  =EFF  (  6  )  *  <  X  <  19  >  ♦X  C  79 )  -MP  <  ',)  -W<  6>  )  *>:  <  1  >  -X  ( 3 1 
G  (  1 9 )  =EFF  (  7 1  v  (  X  i  20 )  «■  X  ( 90  >  +  ( P  (  7 )  -W  ( 7  >  >  *  X  (  i  )  -  X  ( CiO  . 
G  <  20 )  =EFF  <  8 '  M  X  ( 2 1  )  +  X  ( 3 1  )  >  ♦  <  P  <  O )  -  W  <  8  >  )  ►  X  (  I  '  -  X  <  33  ) 
G(21  >=EFF(  9)<  (  X  < 22  )  +  X  ( 82  ))■*•<  P  <  9 )  -H ( 9 ) ) *X  < 1 > -X ( 34 
G(22)»EFF(  10)  *t  X<23>  fX<83)  >  MP<  10>-W<  10>  MX(  1  >-X  ( 
G  <  23  )  =EFF  (  1  1  )*<  X<24>  t-XiS  »  '  1  1  >-W<  1 1  )  )  *X  <  1  >  -  X  ■ 

G  (  24  )  =EFF  (  12 )  •  <  X  ( 25)  >X  (85  >  )  '-IP  (  12  )-W<  12)  >  »x  i  l  )  -X  < 
C 

C  SOIL  DRAINAGE  CAPACITY 

C 

G ( 25  >  =MOri  (  1  >*X<  1  )-X(26> 

G ( 26 ) *MDR ( 2 ) *  X ( 1  •-X(27> 

G(27)=M0R<3)*X<  t )-X(23> 

G ( 23 ) *MDR  <  4 ) «X< 1 )-Xi29) 

G ( 29 ) -NCR ( 3) »x ( 1 ) -X  <  30) 

G( 30) =MDR ( 6) *X ( 1 ) -X<  31) 


.  X  (3). 

X( 13) ) /DP 
i  .  x  (  3 ) , 

X ( 13) ) /DP 
>■ ,  X  t  3 ) , 

X  <  13)  ) /DP 
>,  X  (  3  ) , 
i ( 13) ) /DP 
i,  X(3>, 

X ( 13) ) /DP 
i  ,  X  (  3  ) , 

X< 13) > /DP 
i ,  x  ( 3 ) . 

(•.13)  )  /DP 
>  ,  x  (  3  ) 

=.  •’  13)  )  /DP 
X  (  3  )  , 

••  1  1  3  )  >  'DP 
l  >  x  <  2  ) ,  X  i  3 ) , 
v ; 13))/ DP 
:  <  x  <  p ) ,  x  <  3 )  ■ 
» >’  1  :'■>)/  DP 
;  (  v  <  2  i ,  X  (  3 ) . 
X( 13) ) /DP 


'•  7  > 

j  -t  ) 
37) 
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non 


C  LI  MOD  PROGRAM 


0(31  >='  r:<  ( 7 )  a  < !  n  -  x  ( 33) 

C  (  32  I  =r-cr<  X  (  1  >  -  X  (  33  ) 

G  (33  )  =MDR  .  9  i  *X  (  1  )  -X  (  34  > 
G<31  >=NI \  10»  *  <  (  J  '  -X  (31,) 
G(35)=MLR ( 1 1 >*X ( 1 ) -X(36) 
G(36)=r-U;R  (  13  >*X  (  1  )  -X  (37) 


NITROGEN  BALANCE  AT  IRRIGATION 

31  TL 

G  <  37 )  =  1  0/F 

NL< 1 >*X(62>-< 1 

C-DLN ( 1 ) ) 

»  *(X(  1  4  )  •»  C  X 

•;  :-;o  >  *x  a  2>  +  x  <  49> 

-X (  iZn  > / ( 

A  <  2 )  +  A  (  1 

3'  )  » N I  > 

1  ) *X  (  7 

*  +  <  1.  O-FNL  '■ 

12'  •  /FNL ( 12) A 

"3 )  > +X  <  1  ' 

»X ( 50  ) 

G ( 33  >  =  1  0/F 

NL i 2 ) + X ( 63 ) - ( 1 

0-  L'-CN  (  2  )  ' 

*  *(  X  (  15)  »  (  X  <39)*X  (  3)  t-X<33  ' 

*  v  <  z ) )  /  ( <  ( :< )  +  x  <  2 ) 

) +NI ( 2 , 

.  -X ( 75) 

*  +<1.0  -FNL  < 

1  >  ) /FNL ( 1 ) »<  <>2 

'  )  +  X  (  1  )  *  A 

<  51  ) 

G ( 39 ) = 1 . 0/F 

.ML  <  3  )  *  X  (  64  )  —  i  1 

0  DFN ( 3 ) ) 

*  *  (  X  <  1 6  )  *  (  < 

(  40  )  *X  (  4  )  i-X  <  39 

+  X (3) ) / ( X 

( 4 ) +X  <  3 ) 

t+m  (3: 

.  *  X  (  7  6  ) 

*  +(1  O-FNL'' 

2  >  )  /  FI  JL  <  2  )  «  X  (  63 

>  )  t  X  (  !  )  +  ■ 

A  52) 

G ( 40 )  =  1  0/ F 

N  L  >  4  i  *  X  (  o  5  )  -  >  1 

r-cri.  j  '  ' 

+  *  (  X  (  1  7  )  +  (  a 

>4i ;  *x  ■;  5  >  +x  ;  4  o  ‘ 

A  \  S  ^  .<  /  \  A 

.  3  .  •  X  (  4  ) 

•  *f  N I  %  4  i 

A  (  77  > 

*  +(1  O-FNL > 

3 )  ) /FNL <  3 ) <  X  < 64 

■  )  +  X  (  I  )  *  A 

;  « •  “5  •, 

G  f  41  >  -=1  0/F 

Nl_  .  0  >  *  X  (  6c  )  -  :  1 

DC-  n  -•  3  )  • 

»  *  <  X  <  1  S  >  <M  * 

1 2  )  --  X  <  <■>  *  +  X  >  4 1  ) 

-/.".)./.< 

»'  L  )  >  X  (  5  > 

i  +N I  (  5 

A*(7i.n 

*  +( 1  0-FNl< 

4  )  ■  .-FNL  (4  i  (65 

:•  )  --  x  (  i  >  *  < 

■:  54 

G(42)-l  o. 

•>_  i  6  )  *  a  (  67  i-il 

.  -L-c!  •-•  . 

»  *  (  X  (  I  9  )  *  1  X 

*  >  3  -  *  X i 7  )  ► X ( 42  • 

--  X  (  6  )  •  /  i  A 

(  7  w  X  i  6  ) 

+  i-t  I  (  j 

+  X  ( 7V) 

*  +(1  O-FNL. 

5  J  >  /  -"NL  <  5  )  «  X  >  e  c 

:  )  +  X  l  ;  *  +  » 

.'  55  > 

G  <  43 .  =  1  0/F 

<  v  i  » x  (68 )  -  ( l 

j-C'EN  (  7  >  i 

*  * ( X ( 20 ) *  (  x 

(  44  l  «■  X  »  n  '  ♦‘X  \  43  ' 

+  X  .  7  >  i  /  (  a 

<  a ) + x  <  7  > 

i  f  N I  <  7 

+  X ( 80  > 

*  +( 1  O-FNL ( 

c ) ' /FNL (6) «<(./ 

)  +  X  (  1  1  <•  A 

i  56) 

G  .  4*1 )  -  1  C  -'F 

NL : 3 ' +x <69  '  -  ( l 

‘  r Li.  3  • 

*  *  (  X  (  2  1  >  *  <  V 

i  4 N  >  +  x  (  9  >+/  i  -44 

+  x  :  8  )  '•  -  ';  . 

<  9 ) + x  <  a  > 

-  +M1  ( 

■  -*X  !  81  ) 

*  +(  1  O-FNi.  . 

7  )  ;  /  FNL  (7)  v  X  ;  : 

•  )  f  1  (  :  ;  1  X 

•'  57  ) 

G  <  45  )  =  1 .  C/"7 

fiL  *  ~  *  a  v  7 0  )  -  ••  l 

■  i  -L  ■  !  *  ^  1 

»  - ( X ( 22 )  »a  a 

1  .  10)  N  ;  ■« 

■  •-  X  a  V  ;  i  /  » 

A  l  1  0  )  +  A  . 

4-  ,  )  f-N  I 

-  >  »X  (3. 

*  +(1  O-FNL' 

3)  )  /  FNL  <  3  /  *  X  =•  - 

•  )rX(  1  i  ^  ; 

'  56) 

G( 46 )  =  1  0 / F NL < 10  <  *X . 7  !  •  -  1 

(  1  O’ 

1 

ft  ♦  t 

>  '1 7  i  «•  V  >  1  1  >  -  V  .  1  -> 

‘  A  -.  1  0  )  )  ■ 

.  x  ;  ;  ;  +  •; 

•  0  '  .  - 

A  '  (  J  0  )  + 

*  +<1.0  -FN! 

=  )  ;  •  r  Ni_  (  9  •>  '•  1  • 

'  -  x  (  :  >  -  ■ 

'  ’*  » 

G  (  4  7  )  —  }  .  0 

_  1 1 ) *X ( 72  -  ■  1 

^  •’  “  i  >  1  * 1  ’  1 

fr  ♦  <  X  (  2  4  >  » 

•!6 .  12) -i  X  .  ' 

A  X  1  I  ;  :  ) 

;  i  i  j  .  •  r  » 

A  ' 

'  J  I  »  •=■ 

*  +<1.0  FNL  ■ 

10'  ;  .FNL'  H'.Ni 

"  1  )  )  !  <  >•  1 

»  X  \  6  0  ) 

G ( 43  >  =  1  C  ->F 

Ll  v  1 2  /  •*  X  <  72  •  •  •  * 

O-L'EN  (  12 

)  ) 

*  *  (  X  <  2  5  >  *  <  * 

•  4  7  ;  *■  <  {  13)  i  »  l 

<  X  i  1  2  ;  )  / 

i  i  ,  1  J  )  9*  X 

.  i  k  >  *  * 

ill  12'* 

*  +( 1. 0-FNL< 

1 J  )  > /FNL (11  ’ -  X  > 

7c  /  5  ■♦•X  i  !  > 

+  x  ( 6  J  :• 

CROP  UPTAKE  OF 

un  RCGEN 

G  <  49 )  =  X ( 50 ) 

-CP  <  1  ) »NNU(  1  ! 

*  *  (  1  C  -  E  *  r‘ 

.  (  -NI  ;  1  )  * X  <  74 ) - 

1  0  F NL  • 

12') / FNL 

i  *  i  ■  ^  «*•  / 

.  C.  /  -  %  1 

■2)  > 

*  » <  X  ( 2  >  +a (  1 

:  >  ■  - . .  1 4  > 

»  * ( X ( 38 ) *  X . 

2) +X i 49  >*X (  1 81  ' 

/(  ( X ( 2  )  + 

A  <  1  3  )  ) 

*  * ( X (  l ) +  MNU>  1  I  )  )  ; ) 

G  (  50  )  =X  <  5  1  >-C.P:2'  *MNU(. 


;  (  83  > 

;<£**« 

i  ( B5 ; 
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C  1.1  MOD  PROGRAM 


*  *<  1.  O-EXP  (  ■;  <  -MI  2  )  *X  <  75  >  -  '  1.  0  FNL(  1  >  >/FNL  <  t  XX  it.?)  ) 

*  *(X(3>«-X(2>  )  -X(  15) 

*  »(X(39/»X.3>«-X(25>*X(2)>>''(<X<3)+X(2>) 

»  »(  X  (  1  )  «HMU .  2  )  i  )  )  i 

G  <  5 1  >  =X ( 52 i  -CP  >•  3  )  *MNU(3> 

*  *(  1  O-EXP  <  (  (  -Ml  (  3>  *X  <  V6)-i  1  0-  r-NLi2>  )  /FNL  (2) *X (63) > 

*  *  (  X  (  4  )  -*•  X  i  3  '  )  -  x  (  16' 

*  *  (  X  (  40  )  *X  (  4  )  -rx  <  39  >  *X  (  3  >  )  > /  (  ( X ( 4  >  +X ( 3  > ) 

*  *  (  X  (  i  )  »HNU  '  Oii  '  )  i 

G ( 52 )  =  X ( 53 ' -CP < 4 ) ^MNU<  4  ) 

*  *(  1.  O-EXP  (  (  (  -MI  (4)  *X<77>-(  1  0  -FNL(3>  > /FNL  <  3  )  *X  ( t>4  )  ) 

*  *(XC5)*X(4) )-X( 17) 

*  *<  X(41  >*X  (  5>  «•*  (40)*X  (4)  )  >  /  (  (  X  (  5)+X(4>  ) 

*  *(X(  1  )*MMtJ<4)  )  >  )  > 

G  (  53  )  =X  (  5-t  )  -CP  (  5  )  »NNU  (  5> 

*  »il.  O-EXP  ( ( (-NI ( 5>*X (73>-<  1.  0-  FNL  (  4  )  )  /FNL  (  4  )  *X  ( 65  )  ) 

*  *  (  X  <  6  )  +  X  (  5  )  ) - X ( IS) 

*  *  <  X ( 42  >  *X  <  6 ) ( 4 i ) *X ( 5 i  )  > / ( ( X ( 6 ) +X ( 5 ) ) 

*  *  (  X  <  I  )  *H(  JU  (  5  )  )  )  )  ) 

G  <  54 ) =X  <  55 )  -CP(o)*MNU(6) 

*  *(  1  O-EXP  C  (  (  -NI  (6)  »X  (79)  -(  1  0  FNL  (  5)  )  /FNL  (  5  )  *x  (6o  '•  ) 

*  *  (  X  (  7  )  -t-X  (  6  ) ) —  X  <  19) 

*  *(X (43  MX (7)-*-X ( J2 )*X (o'  '  > / ( ( X (7>  +  X (6) ) 

*  *  ( x  ( i )  *mnu  ( c  > ' )  :•  > 

G ( 55 ) =X ( 56 ) -CP ( 7 ) *NNU ( 7  ) 

*  *  (  1  O-EXP  (  (  (  -M  !  (  7  >  *X  <  BO  >  -  (  1  0— FNL  (  6  )  )  /FNL  (  6  )  •-  X  (  *7  )  ) 

*  *(X(8)+X(7) )-X<20> 

*  *(X(44)*X(S)*X  (43)*X(7  )  )  )  /  (  (  X  (BH-X(7>  ) 

*  *(X( 1 >*KNU(7> )  )  )  ) 

G  (  56  )  =  X  (  57  )  -CP  (  3  >  -  MNiJ  (  8  > 

*  *< 1.  O-EXP (  t  ( — M I ( 3 ) * X  <  6 i  )-( 1.  0~FNL ( 7  > ) /FNL ( 7 > *X ( 68 )  ) 

*  »  (  X  (  9  >  +  X  i  3  '  )  -  X  (  2  1  > 

*  *  (  X  (  4  5 )  *  X  f  9  )  «■ »  (  4  4  )  »  x  ( 1?  '  )  )  /  (  (  X  ( 9 )  +  X  ( 8 )  ) 

»  *  (  X  (  1  )  *nNU  '  Cl '  )  '■  '  ' 

G ( 57 ) =X ( 53  > -CP < 9 >«MNU(9> 

*  *( 1.  O-EXP ( ( ( -MI (9) *X(62>-< 1  0- FNL ( 8  > ) /FNL ( 3  > *X (69  5 > 

*  *(  X  (  10>»X  (  9  !  )  -<  ••  22) 

-  »(X(46'»X  (10 4  0)*X<9))'/((X(I0>  +  X(9)> 

*  •  ( X  (  •  ) »MNU< 9 ) ) '  S 5 

G  i  5  0  '  =  X  (  5  9  ;  -  C  P  (  10)  5  MNU  (  i  0  > 

*  *<  1.  O-EXP  ■,  •  (  -NI  •;  10'  *X  (331  -  <  1 .  0-FNL(9)  )  /FNL  <9)*X(  70  )  ) 

*  *(  X  <  1  1  )  +X  <•  10)  )  “X  (23) 

»  *(  X  (  47>*x  ;  ;  l  >»<  .•  46>*X  (  10)  '  >  /  <  (  X(  1 1  >  *X  <  10)  > 

*  *(  x  <  1  )  V7i>i)) 

G ( 59 )=X  <  60 ’ -CP (11) *MNU( 1 1 ) 

*  ♦(  1.  0-  EXP  >  i  (  -MI  .  1  1  >  *X  <  3-1  ■-(  1  O-rNL  <  10)  )  /FNL  (  1  O'  *X  (  71  )  ) 

*  *( X ( 121+X  »  1 1)  ) -X  <  24 ) 

*  *(  X  <  48>*x  (  i2  >  +  X  (  47)  *X  (  1  li  )  )  /  (  (  X  (  12H-X  (  1  1  >  > 

*  *(X( 1 >*MNU( 11')))) 

G  (  60 )  =X  ( 61  > -CP  (  1  2  )  «  MNW  (  12) 

*  *( 1  O-EXP 1  (  (— NI ( 12) *X (85) -( 1.  O-FNL ( U  > ) /FNL  ( 1  1  ) »X < 7?)  ) 

*  *( X ( 13) +  X ( 12) ) -X ( 25 ) 

*  *<X(49)*X ( 13)+X( 48) *X ( 12)))/((X(13M-X(12) ) 

*  *(X( 1 ) »MNU ( 12) ) ) ) ) 
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C  LI  HOD  PROGRAM 


NITROGEN  BALANCE  IN  STORAGE 

G(61  >  =  X  <  33)  *X  (2) -X  (49  )  *X  <  13  >  «-X  (  14  )  »(  X  <3H)  *X  <2  )  *X  (  49  >*X  (  13)  ) 

»  X ( X ( 2  >  +X (  13)  >♦  5«-NK<  1  >  *<  X  <38>*X  (2>  +  X(4V>  *X  <  13)  >  +NI  (  1  )*X  (74  ) 

G  <  62 >  =  X  (  39  >  *X  (  3  >  —  X  (  33 ) *X  <  2  >  +  X  < 1 5 ) *  <  X  < 39 ) *X ( 3 ) +X  <  38 )  *  X  (  2  )  ) 

*  /  <  X  <3>  *X  v  2  >  >♦  5*NK(2)  »  <  X  v  33)  *X  (8)  *X  (39)  *X  <  3)  )  *N I  ■  2 ) *X  < 75  ) 

C<  A3)  =X  (40>  -X  (4)-X  (39)  *X  <  3  >  *-X  <  16)  *  <  X  (40)  *X  <  4  )  f  X  <39  )  «X  <  3)  ) 

*  /<  X<3)+X<4  > >♦  5»NK(3>*( X  <  40) *X(4  >+X(39>  *X  <3>>  *NI (3) *X (76) 

G  <  64  i  =X  (  4  1  >  *X  (  5  >  -X  (  40  >  *X  (  4  >  <-X  (  1  7  )  ■»  <  X  <  4  1  )  *X  (  5  >  +X  (  40  ;  <  X  ( <’,  )  ) 

*  /  <  X  <  5  >*-X  <  4  .  )  f.  5»NK <  4  >»  <  X  .  40)  *X  <  4>  +  X  <41  )  *X  <  5)  >  «-NI  <  4  )  *X  <  77  ) 
G<65X=X<42  )  »X  <6)-X  <41  )*X  (  5  )*X  <  18)  *<  X  <42)*X  <65  *X  <  41  )*X  <  5)  ) 

*  /<X<6>t-X(5>)+  5»NK<5>*  <  X  •  41  )  *X  <  5>+X(42)  *X  <6)  )  »NI  <  5)  *X  (70) 
G(66)«X(43)»X(  7>-X  (42>*X  <  6)  <-X(  19>*(  X(43)*X(7)  <-X(  42)*X  (6)  > 

*  /  (X(7)+X<6)  >•*•  5*NK(o)  *(  X  <  43  >  *X  (  7  WX  <  42  >  *X  <  6  >  > +N I  ( 6  >  *X  <  7? ) 

G(67)  =  X (44  . <  8)-X (43) *X ( 7 ' • X  <20)*<  X (44) *X  <S ) <  X  <  43 i *X  <  7 ) > 

*  /  <  X  (8)  +X<  7'.  )  +  5*NK<7)  *(  X  <  44)  *X(8)+X(43)*X  (7)  )-rNI  »  7)  *X  (SO) 
G<60)=>X  (4>5  ,-»X  <  9)  -  X  (44)  *X  (3)  i-X  <21  )*<  X  <45)*X  (9)  frX  <44  )*X  <t!)  ) 

*  /  (  X<9)  »X(3)  )  *  5*NK(S)*(  X  ••  45)  ifX  <9)+X  (44  )-*X<e>  )+NI<  9)  *X  <ei  ) 
G(69)=X(4t ,  s-x  <  10)-X<45)  »X  > .9)  tX(22)*<  X  <  46)  »X  <  10J+X  <  45  )  *X  <9  )  ) 

*  /(X<10)*X<9))*  5-s-NK  (  9  )  *  <  x  >  45  )  <  X  (  9  )  -s-X  ( 46 )  4  X  i  1 0  )  )  *N  I  <  9  )<•  X  <  U 
G(70)=X(47) -  X<  i  1  )-X(46)  *X  :  10)  <-X(23)*(X<  17)*:<  :  1  1  )+x  ,  46)  -  Xt  J 

*  /<X<11)+XU0)'>+  5*NK(  10)  -><  X  (47)  4-X  <  11  )  ’■X  <46>»X  <  10  >  )  +NI  <  10)  *X  <  S3) 
G  (  7 1  )  =X  (  43  >  *  X  (  1 2  >  -X  <  47  )  *  x  :  1  1  )  +X  ( 24  )  *  (  X  (  43  )  «  X  <  J  2  )  '►x  <  47  )  «■  X  <  1  1  )  ) 

*  /<  X  <  12)  «■'<  <  1 1  )  >+.  5»NK  (  11  ;■  *<  X  <43)*X(  12)  +X  <47  >  ■=  X  i  1  1  ;  1  N I  <  1  1  )  *?,  (  84  ) 

G(72)=X  <49  >  *X  (  13>-X  (  40)  *  <  •  1  2  >  --X  <  25  )  *  <  X  <  49 )  *  <  1  3  )  +  '  <  4B  )  *  X  (  )  2  1  ) 

*  /  (X  <  12)-»-X  (  13)  >♦  5*NK(  12*  *  <  X  <4B>  *X<  12>fX  (49)  »X  <  13)  >*IMI  <  12)  *X  <  35) 

NITROGEN  CONC  IN  PERCOLATE  CONSTRAINT 

G  <  73  )  =NCQ*  <  X  <  26 ) -X  <  27  )  t-X  <  2S )  +X  (2? ) 

*•  +  X  ( 30  >  *X  (  31  )  *-X  <  32  >  «-X  ( 33  > -X ( 34  > 

*  +X  <35)  +  X  (36M-X  (37)  )  —  <  X  (62  H  X  (63)+X(64) 

*  +X ( 65 ) +  X ( 66 >-T < 67 ) +X ( 6S ) *  X ( 6? ) 

*  +X  (70  )  +X  (  71  )  *-X  (  72  >  +X  (  73  )  ) 

MONTHLY  QUALITY  OF  PERCOLATE  CONSTRAINTS 


G ( 74 ) =NCGM ( 1 >*X<26) -X<62) 

G  (  75  )  -WCCM  (  2  >  *X  (  2  7  )  -X  ( 63  > 
G(76)=NCuM(3)*X(29)-X(64 ) 

G ( 77 ) =NCGM  <  4 ) * x (  29 ) -X ( 65  > 

G ( 78 ) =NCQM <  5  >  * X  <  30  > -X ( 66  > 

G  (  79  )  =NCGM  <  6  )  *X  (  31  ) -X  (  67  '• 

G ( eo ) =NCGM <  7 ) 1 1 ( 32 ) -X  <  68 i 
G  (31  )  =NCGM  i  3  )  -»X  <  33  )  -X  (  69  ) 

G ( 32 ) *NCGM ! 9 ) *X  <  3  4 )-X ( 70  > 

G ( 33 ) =NCQf  1  ( 1C)*X(35>— X(71  ) 
G(34)=NCGM( 1 1 ) *X  <  36)  -X  <  72 ; 
G ( 35 ) =NCGM  < 1 2 ) »X ( 37 ) -X ( 73 ) 

OBJECTIVE  FUNCTION 


G(36)=0A*X ( 1 ) «-OB*DMAXl ( X  <2  - , X(3). X(4), X(5>. X(6> , 


O  PJ 
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LI  MOD  PROGRAM 


*  X  (  7  > .  X(S>.  X  l  9  ) ,  X  i  10).  X  (  l  1  >  •  X ( 1 2 ) ) 

NITROGEN  REMOVED  IN  STORAGE 

C(37>=  5*NK< 1 )*<X(08)*X<2)+X(4V)*X<13) ) 

*  +.  5*NK<2)  » (  X  ( 39  >  »  X ( 3 ) +  X  <  33 ) *  X ( 2 )  ) 

*  +  .  5*NK<3>  »(  X  (40)  *X  <  4  )  *'<  <  39  )  *X  (3)  ) 

*  +.  5*NK<  4  >  *  <  X  (  4  1  >-*X<5)+X<40)*X<4)  > 

*  +  5*NK<  5>-*<  X(42)»X<6)+X(41  >  *  X  (  5  )  ) 

*  5*NK<6>*<  X  (  40  W<  (7  )  +X  (•  4  4  )  *  X  (6)  ) 

*  +  .  5*NK<7>  *<  X(44';-»X(8)»X(43)*X(7)  ) 

*  +. 5*NK<Q)*iX(45)»X<9>+X<44)*X<8) ) 

*  +. 5*NK ( 9  >  ♦  < X (46)*X ( 10)+X (45) *X (9) ) 

*  +.  5*NK< 10) * ( X (47) *X  < 1 1  )  *XM6)*X(  10)  > 

*  +.  5*NK(  1  1  >  »<  X  <43)*X<  12) +X  >'47)«X(  1  1  i  ) 

*  +.  5*NK< i2>*<  X (49) »X< 13) +X (48)*X ( 12) ) 

PERCENTAGE  OP  NITROGEN  REMOVED  IN  STORAGE 

G(88)=G<87>  /  <NI  i  1  )  *1  (  1  )+NI  (2)  *1  <2)+NI  (3)*I  (3)+NI  (  4  >  «•  I  (4) 

*  +N I  (  5 )  *  I  <  5  <  +N I  ( 6  )  *■  I  ( 6  )  +N I  (  7  i  *  I  (  7  )  +N I  ( 8 )  *  I  (  8  ) 

*  +NI (9 >*I ( 9 ) +N I < 10) *1 < 10)»NI (11 ) *I< 11 )+NI ( 12>*I ( 12) ) 

NITROGEN  REMOVED  BY  CP0P<S) 

G(09)  =  X<  l  )  *v  X<  50)fX<  Sl)+X<  32)  «-X(53)  +X(54)+X  (5b) 

*  +  X ( 56 ) +X ( 57  >  +  X  <  58 ) +X  <  57  >  ’■X  ( SO )  +X  <61  )  ) 

PERCENTAGE  OP  NITROGEN  REMOVED  BY  CRDP(S) 

G(90)=G<89>  /  <NI  (  1  )*I  (  1  >+NI  (2)  *I  (2)+NI  (3)*I  <3>«NI  <  4X1  <4> 

*  +NI  <5>*I  <  5  '-"MI  <6>*I  <6)i-NI  <7>*I  <7>+NI  (B)*I  <8> 

*  +NI  (9)*I  (9)-!NI  ( I0)*I  (  lO’-Nl  U  1  )  *1  ( 1 1  )+NI  ( 12i*I  <  12)  ) 

NITROGEN  REMOVED  IN  SYSTEM 
G  <  9 1 ) =G  <  37 ) +G  <  89 ) 

PERCENT  QF  NITROGEN  REMOVED  3Y  SYSTEM 

G<92)»G(91  )  /  <NI  <  1  >«  I  (  1  )  f-NI  <2)  »I  (2)+NI  (3)*I  <3>*NI  (4)*I  (4  ) 

*  +NI<5)*I<5)+NI(6)eI(6)+NI < 7 > *  I < 7 > +NI < 8 > *1 ' ~ > 

*  +NI ( 9 ) *  I <  9 ) +NI < 10) *1 < 10)*NI < 11 ) *1 ( 1 1 )+NI ( 12)*I< 12) ) 

MAXIMUM  WEEKLY  APPLICATION  RATE 

G  <  93 ) =DMAX 1 ( X ( 1 4 ) >  X  (  1  5  ) ,  X  < 1 6 ) •  X ( 1 7 ) .  X ( 10).  X ( 1 7 ) , 

*  X<20) , X<21 )+X(22)+X(23>+X(24)+X (25) >/<X( 1 )*4.  3) 

END 
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